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As  systems  manager  for  the  Urban  Mass  Transportation  Adminis- 
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-Transportation  Systems  Center  (TSC)  is  conducting  research  and 
development  efforts  directed  toward  the  introduction  of  improved 
technology  in  urban  rail  system  applications.  As  part  of  this 
“program,  TSC  is  conducting  analytical  and  experimental  studies 
toward  improved  safety -in  urban  rail  systems.  A specific  goal  in 
this  area  of  safety  is  to  reduce  the  number  of  injuries  that  may 
result  from  the  collision  of  two  trains. 
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four  volumes.  Part  1 describes  the  results  of  Task  1 which  is 
concerned  with  the  initial  impact  of  two  transit  cars.  The  results 
of  Task  2 which  is  concerned  with  the  primary  collision  of  two 
impacting  transit  car  consists  are  described  in  Part  2.  Part  3 
describes  the  results  of  Tasks  3 and  4 of  this  study  which  are 
concerned  with  prediction  of  passenger  injury  and  guidelines  for 
evaluation  of  railcar  designs.  The  final  volume  is  a manual 
containing  a description  of  the  organization  and  use  of  the  IITRAIN 
computer  code  which  was  developed  as  a tool  to  help  meet  the  goals 
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1.  INTRODUCTION 


The  collision  of  two  consists  of  transit  cars  ca_n  be  broken 
into  three  separate,  but  interdependent,  phenomena:  initial  impact, 
primary  collision,  and  secondary  collision.  Initial  impact  is 
concerned  with  the  mechanics  of  the  initial  impact  of  the  leading 
cars  of  two  consists.  The  interaction  of  all  of  the  cars  and  car 
components  of  two  impacting  consists  comprise  the  primary  colli- 
sion. Secondary  collisions  include  the  interaction  of  passengers 
with  the  car  components,  passengers  with  passengers  and  passengers 
with  other  loose  objects.  This  final  report,  submitted  in  four 
volumes,  describes  the  results  of  the  IIT  Research  Institute 
(IITRI)  program  which  is  concerned  with  the  collision  of  transit 
car  consists  on  straight  level  track.  Part  1 of  the  final  report  is 
concerned  with  the  initial  impact  of  the  leading  cars  of  two 
consists.  The  results  of  the  study  of  the  primary  collision  of  two 
impacting  consists  are  given  in  Part  2,  and  Part  3 is  concerned  with 
secondary  collisions  including  the  prediction  of  passenger  injury 
and  guidelines  for  evaluation  of  new  railcar  designs.  The  final 
volume  is  a manual  containing  a description  of  the  organization  and 
use  of  the  IITRAIN  computer  code  which  was  developed  as  a tool  to 
help  meet  the  goals  of  this  contract. 

1.1  Program  Objectives 

The  program  objectives,  as  taken  from  the  contract,  are 
restated  here. 

Item  la:  Formulate  an  analytical  model  in  two  dimensions, 

longitudinal  and  vertical,  of  the  leading  cars  of  two  impacting 
consists  in  sufficient  detail  to  examine  the  mechanics  of  head-on 
initial  impact  on  straight  track.  This  model  will  include  the 
distribution  of  mass  in  the  cars  as  well  as  the  nonlinear  force- 
deformation  relationships  existing  among  major  structural  subas- 
semblages. Consideration  will  be  given  to  the  shapes  and  config- 
urations of  the  impacting  surfaces  and  to  the  forces  generated  by 
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the  impact.  The  model  shall  be  capable  of  establishing  the 
critical  parameters  which  govern  whether  the  cars  crush,  displace 
vertically  and  override,  or  crush  with  subsequent  override. 

Item  lb:  Utilize  the  above  analytical  model  of  initial  impact 
to  assess  impact  controlling  devices  currently  in  service,  such  as 
-anti climbers , couplers  and  draft  gears  of  various  designs.  This 
assessment  shall  uncover  the  critical  parameters  of.  such  devices 
which  govern  whether  the  cars  crush,  displace  vertically  and  over- 
ride or  crush  with  subsequent  override.  The  contractor  shall 
develop  recommendations  concerning  future  directions  of  effort  in 
design  of  impact  controlling  devices  which  would  be  particularly 
pertinent  to  crashworthiness  goals. 

Item  lc:  Develop  an  experimental  test  plan  for  the  evaluation 
of  the  strength  and  effectiveness  of  future  designs  for  impact 
controlling  devices.  These  tests  are  to  assure  that  the  forces 

generated  during  impact  do  not  produce  structural  failure  of  the 
impact  controlling  device  or  vertical  misalignment  and  override  of 
the  car  body.  The  test  plan  is  to  be  sufficiently  detailed  so  that 
all  equipment,  fixtures,  instrumentation  and  procedures  are 
completely  described. 

Item  2a:  Develop  an  analytical  model  in  two  dimensions,  longi- 
tudinal and  vertical,  of  the  primary  collision  of  two  impacting 
consists  of  urban  railcars  of  similar  and  different  configurations. 
This  model  will  include  the  formulation  of  the  leading  cars  devel- 
oped in  Part  1 of  this  program,  as  well  as  the  distributions  of  mass 
and  nonlinear  force-deformation  relationships  existing  among 
major  structural  subassemblages.  This  model  shall  be  capable  of 
determining  the  extent  of  crushing  and/or  override  suffered  by  the 
individual  cars  in  the  consists,  as  well  as  the  time  histories  of 
displacement,  velocity,  and  acceleration  in  both  the  longitudinal 
and  vertical  directions. 

Item  2b:  Develop  methods  for  generating  the  dynamic  force- 

deformation  relationships  for  structural  subassemblages  comprising 
the  critical  modules  of  railcars.  These  methods  shall  include 
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finite-element  analysis,  scale  modeling  and  full-scale  testing 
procedures  including  specifications  for  required  testing  equipment 
and  instrumentation.  Utilize  the  finite-element  analytical  method 
to  generate  the  nonlinear  force-deformation  relationships  among 
major  components  of  a typical  urban  railcar. 

Item  3:  Develop  the  analytical  methodology  of  passenger 

injury  due  to  secondary  collision  to  include  modes  of  injury  due  to 
longitudinal,  vertical,  and  pitching  motions  of  the  vehicles  after 
impact.  This  methodology  shall  be  capable  of  considering  the  loca- 
tion of  the  passenger  prior  to  impact,  his  orientation  (seated, 
standing,  facing  forward,  facing  sideways,  facing  rearward),  the 
configuration  of  interior  features  of  the  cars,  passengers  density, 
and  passenger  restraint.  This  methodology  shall  also  be  capable  of 
determining  the  severity  of  the  injury  sustained  by  the  passenger. 

Item  4:  Utilize  the  results  of  Items  1 through  3 to  develop 

guidelines  for  the  evaluation  of  proposed  railcar  designs,  and 
guidelines  for  the  development  of  new  railcars.  These  guidelines 
are  to  be  developed  in  parametric  form,  so  that  individual  param- 
eters may  be  considered  and  the  effects  of  specific  values  assigned 
or  computed  for  these  parameters  may  be  assessed.  These  parameters 
are  to  include: 

a - the  number  of  cars  in  the  consist 

b - operational  velocity  ranges 

c - dimensions  and  weights  of  each  car 

d - placement  and  dimensions  of  windows  and  doors 

e - placement  and  weights  of  mechanical/electrical 
equipment 

f - interior  configurations  of  passenger  compartment 

g - carbody  force-deformation  relationships  between 
major  structural  subassemblages 

h - locations  of  carbody  centers  of  gravity  (c.g.) 
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1.2  Passenger  Survivability  Considerations 

The  objective  of  the  task  reported  in  this  volume  was  to 
assess  the  nature  and  severity  of  potential  passenger  injury  and 
fatality  producing  mechanisms  due  to  secondary  collisions  which 
occur  in  intracity  train  crash  environments  and  to  use  this  analyt- 
ical methodology  to  identify  the  factors  and  car  design  parameters , 
including  car  interior  components  which  significantly  influence 
these  events.  In  this  manner,  meaningful  design  guidelines  can  be 
established  which  will  minimize  and/or  limit  the  hazard  to  which 
the  passengers  will  be  exposed  in  the  event  of  most  crash  environ- 
ments . 

The  overall  problem  is  very  complex;  however,  a rational  exam- 
ination of  the  many  aspects  of  the  phenomena  involved  should 
provide  an  insight  into  the  establishment  of  practical  steps  which 
can  be  taken  to  mitigate  the  most  undesirable  events.  The  train 
crash  conditions  of  interest  may  also  be  varied,  but  clearly  the 
major  change  in  motion  is  colinear,  corresponding  to  main  frame 
decelerations  (or  accelerations  for  certain  cases  in  a consist- 
consist  collision)  resulting  in  a specified  change  in  gross 
velocity.  Since  only  secondary  collisions  are  of  interest  for  this 
task  any  "crushing"  environments  involving  passengers  are  not 
considered.  The  injury  producing  mechanisms  of  concern  are  those 
resulting  from  the  relative  motion  of  the  car  or  car  component  and 
the  passenger  or  passengers  and  their  subsequent  interactions.  The 
interior  configurations  of  the  cars  which  comprise  the  seats,  seat 
arrangements,  stanchions,  etc.,  represent  rather  complex  and 
sometimes  relatively  rigid  boundaries  or  surfaces  at  which  many 
severe  interactions  will  occur.  The  motion  of  the  passenger  will 
be  as  varied  as  the  range  of  passenger  positions,  orientations, 
passenger  alertness,  etc.  The  details  of  a passenger-structure 
interaction,  including  the  structural  response  of  the  local  car 
component  and  the  biomechanical  and  related  physiological  response 
of  the  body  or  body  component,  will  define  the  nature  and  severity 
of  any  secondary  collision  induced  injury.  Thus  the  approach 
selected  was  to  establish  injury  and/or  fatality  criteria  for  a 
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variety  of  body  component  categories,  to  examine  the  response 
characteristics  for  a class  of  both  simple  and  complex  structural 
systems  common  to  train  cars,  and  to  determine  the  motion  details 
of  both  passenger  and  car  components  for  the  spectrum  of  crash 
conditions  of  general  interest.  These  three  basic  elements  were 
assembled  into  a computer  code  with  which  to  evaluate  injury 
profiles  for  specific  train  crash  environments  (established  from 
consist  crash  analysis  such  as  IITRAIN),  passenger  distributions, 
passenger  density,  interior  car  component  properties,  and  interior 
configurations.  Statistical  anthropometric  data  were  used  along 
with  other  uncertainty  factors  to  yield  results  which  reflect  real 
world  variances  which  are  expected  to  be  present  in  an  actual  crash 
situation. 

The  multiplicity  of  situations  which  can  occur  is  so  great 

that  if  one  is  to  proceed  in  an  effective  manner  it  is  essential  to 

identify  the  more  common  and  more  hazardous  ones  and  to  concentrate 

the  available  resources  in  these  areas.  An  onboard  accident 

scenario  for  intracity  train  crashes  was  not  established;  however, 

some  insight  into  such  a scenario  may  be  obtained  by  examining  a 

1 

somewhat  related  scenario  for  transit  bus  operations.  Such  a 
scenario  is  presented  in  Figure  1.  The  bus  motions  which  produce, 
by  far,  the  most  onboard  accidents  are  those  relating  to  decelera- 
tions and  accelerations  and  thus  correspond  to  the  general  type  of 
vehicle  motion  expected  from  consist-consist  or  consist-abutment 
collisions.  To  a certain  extent  this  similarity  justifies  the  use 
of  the  entire  bus  scenario  for  the  intracity  train  crash  situation. 
Clearly  the  most  vulnerable  passengers  are  those  which  are  standing 
in  a relatively  open  area.  Furthermore,  the  standing  passenger 
very  likely  will  not  make  use  of  any  assist  devices  and  an  injury, 
if  it  does  occur,  is  likely  to  be  the  result  of  falling  to  the  floor 
or  the  train  car.  Note  that  the  categories  "female”,  age  group 
"over  50",  and  especially  "over  65",  are  by  far  the  most  common, 


Mateyka,  J.  A.,  "Safety  Considerations  in  Design  of  New  Transit 
Bus  Seats",  Proceedings  of  the  Eighteenth  Stapp  Car  Crash 
Conference,  p 71,  Society  of  Automotive  Engineers,  1974. 
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BUS 

MOTION  AT  TIME  OF  ACCIDENT 

PASSENGER  LOCATION 

. • 

56$  Decelerating! 

• 

39$  Forward  of  First  Cross  Seat 

• 

21$  Normal  Operation 

• 

32$  First  Cross  Seat  to  Rear  Door 

“ ‘ • 

• 

16$  Accelerating 
7%  Turning 

• 

25$  Behind  Rear  Door 

PASSENGER  AT  TIME  OF  ACCIDENT 

WAS 

PASSENGER  CARRYING  OBJECT? 

• 

46$  Standing 

• 

|54$  Yes 

• 

• 

30 % Sitting 
17%  Walking 

• 

46$  No 

• 

7%  Unknown 

WHA1 

• 

WAS  OBJECT  BEING  CARRIED? 
47$  Package 

PASSENGER  USE  OF  ASSIST  DEVICES 

• 

33$  Purse 

• 

28 % Yes 

• 

14$  Umbrella 

• 

|72%  No 

• 

6$  Child 

WHICH  DEVICE  USED?  . 

SEX 

OF  INJURED  PASSENGER 

• 

34$  Stanchion 

• 

[82$  Female 

• 

• 

15 1 % Seat  Handle 
5%  Overhead  Bar 

• 

18$  Male 

HOW 

INJURED? 

AGE 

GROUP  OF  INJURED  PASSENGER 

• 

!6 1 % Fell  to  Floor 

• 

153%  Over  50 

• 

17$  Hit  Seat 

• 

47$  Under  50 

• 

• 

12$  Hit  Stanchion 
9$  Hit  Farebox 

• 

18$  Over  65 

• 3$  Hit  Driver  Partition 


FIGURE  1.  ONBOARD  ACCIDENT  SCENARIO  FOR  CURRENT 
BUS  OPERATIONS  From  Mateyka  (Ref.  l) 
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suggesting  that  physical  alertness  or  responsiveness  are  important 
factors  in  the  occurrence  of  injury  or  fatalities  in  typical  train 
crash  environments. 

1.3  Passenger  Configurations 

Clearly  the  standing  passenger  should  be  given  special  atten- 
tion in  any  meaningful  methodology  development  effort.  In  addi- 
tion, the  passenger  configurations  of  forward  facing  seated 
passenger,  backward  facing  seated  passenger,  and  sideward  facing 
seated  passenger  are  of  great  interest  and  have  been  included  in 
the  current  methodology.  These  four  basic  passenger  configurations 
are  illustrated  in  Figure  2.  The  difference  between  the  forward 
facing  and  the  backward  facing  seated  passenger  is  the  direction  in 
which  the  precrash  velocity  changes.  The  forward  facing  seated 
passengers  will  be  thrown  forward  tending  to  leave  the  seats 
whereas  the  backward  facing  seated  passengers  will  be  pressed  back 
into  their  seats. 
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(c)  Sideward  Facing  Seated  Passenger 
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FIGURE  2.  PASSENGER  CONFIGURATIONS 


2.  BIOMECHANICAL  RESPONSE  TO  IMPACT 


2.1  Basic  Considerations 

Data  for  injuries  resulting  from  falls,  automotive  accidents, 

and  other  related  secondary  collision  accidents  have  been  collected 

2 

and  analyzed  for  many  years.  In  addition,  human  tolerance  experi- 
ments to  various  motion  excitations  have  been  conducted.  Thus,  a 
vast  body  of  data  exists  upon  which  to  estimate  the  nature  and 
degree  of  an  injury  produced  by  a given  impact  stimulus.  Unfortu- 
nately, in  the  case  of  actual  injuries,  the  preconditions  of  the 
events,  such  as  the  impact  velocity,  are  frequently  poorly  defined. 


It  is  clear  that  the  physical  condition  of  the  individual 
involved  in  an  accident  as  influenced  by  age,  diet,  etc.,  has  a 
significant  relationship  upon  the  details  and  consequences  of  the 
event.  Thus  data  acquired  with  the  use  of  cadavers  are  in  many 
instances  biased  in  some  manner.  This  general  data  base  has  been 
examined  where  relevant  and  used  in  the  subsequent  methodology. 


A number  of  injury  criteria,  actual  injury-fatality  criteria, 

have  been  established  which  can  be  used  to  define  the  nature  and 

degree  of  injury  to  whole  body,  head  or  other  body  components  due  to 

impact  with  blunt  objects  or  surfaces,  i.e.,  under  nonpenetrating 

conditions.  While  these  criteria  are  relatively  crude  and  somewhat 
* 

uncertain,  they  represent  the  current  state  of  the  art  in  injury 
criteria  and  are  probably  adequate  for  the  initial  development  of 
an  injury  prediction  methodology. 


The  most  sophisticated  and  documented  impact  severity  criteria 
is  the  Gadd  Severity  Index  (SI)  or  variants  thereof.  This  index  is 
defined  as 


SI  = 


[a  ( t ) ] n dt 


(1) 


Cassidy,  R.  J.  and  Romeo,  D.  J. , "Assessment  of  Crashworthiness  of 
Existing  Urban  Rail  Vehicles",  Calspan  Corporation  Final  Report 
ZL-5368-K-1,  June  1975. 
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where 


t = time  (in  sec) 

a(t)  = the  acceleration  history  of  the  body  or  body  component 
(measured  in  g units) 

n = an  exponent  which  depends  upon  body  component 

t^  = duration  of  the  acceleration  pulse, 

Hote  that  the  change  in  velocity  of  the  body  or  body  component,  VQ, 
due  to  the  acceleration  pulse  is 

ftd 

VQ  =J  a ( t ) dt  (2) 

0 

Thus  when  n=1  and  there  is  no  rebound,  the  severity  index  is  propor- 
tional to  the  incident  velocity.  The  acceleration  pulse  a(t)  is  a 
relatively  strong  vector  function  of  position  on  the  body  (compo- 
nent) and  thus  a specific  location  should  be  considered  to  make  the 
above  definition  of  a severity  index  meaningful.  Such  refinements 

have  not  been  accomplished  to  date!  The  exponent  n has  been  evalu- 
-3 

ated-'  for  frontal  head  impact  conditions  and  a value  of  2.5  best 
fits  the  limited  set  of  data.  This  value  suggests  the  increased 
importance  of  the  peak  or  the  near  peak  value  of  the  acceleration  in 
contributing  to  the  severity  of  a head  impact  event. 

It  has  also  been  determined  that  the  duration  of  the  accelera- 
tion pulse  may  define,  in  addition  to  the  value  of  the  severity 
index,  the  nature  of  the  injury,  such  as  skull  fracture  versus  a 
concussion  for  some  head  impact  situations.  Peak  accelerations 
have  also  been  correlated  to  some  injury  data.  The  value  of  n for. 
use  in  evaluating  the  severity  index  for  other  body  components  has 
been  attempted  but  these  are  estimates  and  a few  firm  data  are 
available  with  which  to  support  the  estimates. 

The  response  of  some  body  components  to  severe  motion  stimuli, 
such  as  the  thorax  and  neck,  have  been  made  but  they  are  of  rather 
limited  value.  The  data  for  the  flexion  (dorsal)  response  of  the 
neck  were  used  in  a limited  manner  and  proved  to  be  quite  useful 
(see  Section  4.7). 

_ 

Snyder,  R.  G.,  "Human  Impact  Tolerance",  Highway  Safety  Research 
Institute,  The  University  of  Michigan,  SAE  Preprint  700398,  1970. 
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2.2  Injury-Fatality  Criteria 


The  goal  of  the  injury  prediction  methodology  is  to  predict 
the  number,  nature,  and  severity  of  injuries  which  occur  in  a given 
intracity  train  crash  situation.  The  severity  of  an  injury  can  be 
expressed  in  terms  of  an  injury  scale  ranging  from  "minor"  to 
"critical"  and  being  bounded  externally  by  the  categories  "no 
injury"  and  "fatal".  These  categories  are  listed  in  Figure  3* 

Their  definitions,  at  least  for  single  injuries,  are  generally 
accepted  within  the  medical  community.  The  precision  with  which  a 
given  severity  index  can  be  related  for  this  injury  scale  is  gener- 
ally poor.  It  appears  reasonable  to  use  this  many  levels  on  the 

injury  scale  as  a means  of  obtaining  a more  detailed  description  of 

the  outcome  of  a postulated  train  crash  such  that  at  least  the 

relative  influence  of  certain  car  design  features  can  be  evaluated. 
The  problem  then  is  to  select  one  or  more  injury-fatality  criteria 
relationships  applicable  to  a given  impact/body  component  mode  and 
to  relate  it  to  the  injury  scale.  This  step  has  been  accomplished 
and  these  relationships  are  presented  in  Figure  3-  They  have  been 
used  in  the  injury  prediction  methodology.  It  is  clear  that  these 
specifics  are  not  well  founded  in  facts;  nonetheless,  they  are  the 
best  available  (or  at  least,  they  are  what  is  needed).  Hopefully 
these  specifics  are  reasonably  accurate  in  a few  instances,  and  an 
attempt  was  made  to  make  them  relatively  consistent.  When  better 
data  become  available  these  criteria  can  be  updated. 

A severity  index  value  of  1000  sec  for  a frontal  head  impact 
generally  corresponds  to  an  injury  in  the  minor-to-moderate  range 
whereas  a value  of  2000  sec  suggests  a cr i t i cal-to-f atal  injury 
level.  These  two  points  were  used  to  establish  the  criteria  curve 
for  a head  impact  and  related  to  a Gadd  Severity  Index  for  the 
condition  of  n=2.5.  Since  the  value  of  n for  the  other  body  compo- 
nents or  impact  locations  such  as  the  face  are  not  well  defined,  a 
value  of  n = 2.5  was  selected.  In  this  manner,  the  ' calculation 
procedure  used  in  the  methodology  would  be  simplified.  It  would 
appear  that  the  general  trend  implicit  in  the  Gadd  Severity  Index 
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FIGURE  3.  INJURY  SCALE  FOR  IMPACT  OF  BODY  COMPONENTS 


form  is  reasonable  and  that  n would  certainly  be  greater  than  unity 
(a  simple  velocity  dependence).  Other  body  components,  such  as  the 
chest,  are  most  likely  not  quite  as  sensitive  to  acceleration  peaks 
as  the  head  would  be,  hence  the  value  of  n is  probably  somewhat  less 
than  2.5.  Nonetheless,  such  a value  was  used  throughout.  The 
criteria  for  the  face  were  selected  on  the  basis  that  at  low  levels 
-of  excitation,  more  damage  would  occur  to  the  face  than  to  the  head. 
The  opposite  behavior  was  believed  to  be  true  at  high  levels  of 
excitation.  The  face  structure  is  thus  envisioned  as  an  energy 
absorbing  system  which  starts  to  crush  at  low  stress  levels  and 
isolates  the  sensitive  brain  at  higher  stress  levels.  The  curve 
for  the  chest  was  placed  well  below  that  for  the  face  in  view  of  its 
(apparent)  stronger  structural  characteristics. 

It  should  be  obvious  that  the  injury  criteria  currently 
available  greatly  limit  the  accuracy  that  any  injury  prediction 
methodology  can  make.  This  aspect  of  the  overall,  rather  complex 
problem  is  perhaps  the  most  limiting  factor  encountered  in  the 
development  of  the  methodology. 
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3.  IMPACTED  SURFACE/STRUCTURE  EFFECTIVENESS 


3.1  Basic  Considerations 

It  has  been  established  that  the  level  of  injury  vthat  a body 
component  sustains  upon  impact  with  a relatively  rigid  object  can 
be  characterized  by  some  aspect  of  the  "lumped  mass”  acceleration 
details  experienced  by  the  body  component  or  parts  thereof  during 
the  impact  process.  It  is  clear  that  the  acceleration  details  and 
therefore  the  injury  level  will  be  influenced,  perhaps  significant- 
ly, by  the  nature  and  responsiveness  of  the  object  (i.e.,  the 
surface  and/or  structure),  especially  at  the  contact  point.  The 
response  of  car  interior  components,  such  as  a seatback,  will 
depend  upon  the  totality  of  imposed  loads,  their  points  of  applica- 
tion,  their  time  phasing,  and  the  displacement  histories  of  the 
attachment  points  of  the  structural  component.  Thus  the  influence 
of  the  response  characteristic  of  the  impacted  surface/structure 
upon  the  injury  level  may  be  quite  complex. 

A series  of  mechanical  response  functions  was  selected  to 
represent  simple  yet  typical  impact/surface  responses.,  and  the 
acceleration  histories  for  a simple  lumped  mass  impact  situation 
were  evaluated.  These  acceleration  histories  were  then  used  to 
determine  the  corresponding  injury-criteria  (index)  dependence  and 
sensitivity  on  the  structural  and  impact  parameter  values.  The 
sensitivity  analysis  is  presented  in  terms  of  influence  coeffi- 
cients and  demonstrates  that  the  response  function  examined,  the 
incident  or  impact  velocity  of  the  passenger  or  body  component  is  a 
relatively  influential  variable.  It  suggests  that  significant 
reduction  in  injury  can  be  achieved  by  retarding  the  relative 
motion  between  the  passenger  and  the  train  component  prior  to  the 
local  impact.  A brief  examination  of  a two-component  structural 
system  was  performed  to  illustrate  the  significant  role  that 
attachment  stiffness  can  play  in  the  design  of  car  interior  sub- 
structures or  components. 
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A potential  load-deflection  relationship  is  illustrated  in 
Figure  4 for  the  local  deformation  at  the  impact  point  of  a body 
component  of  weight,  W . The  impact  is  assumed  to  be  normal  to  the 
surface  and  the  relative  incident  or  impact  velocity  is  V . The 
soft  surface  layer  of  the  body  component  and/or  the  existence  of  a 
thin  layer  of  padding  on  the  surface  will,  in  this  analysis, 
■provide  for  a somewhat  distributed  contact  or  application  area,  A . 
The  initial  resistance  after  the  body  component  and  surface  are 
effectively  mated  will,  most  likely,  be  a relatively  stiff  and 
elastic  resistance.  If  sufficient  energy  (i.e.,  kinetic  energy  of 
the  body  component)  is  available  then  the  structure  may  yield  in 
some  fashion,  more  probably  with  some  strain  hardening  effect  until 
all  of  the  incident  energy  is  absorbed  by  the  surface/structure. 
At  this  point  the  body  component  motion  will  be  fully  arrested  and 
some  bounce  or  throwback  effect  will  subsequently  occur  as  the  load 
drops  during  the  elastic  unloading  phase.  Finally,  the  pad  or 
surface  layer  will  unload.  The  area  under  the  load  deflection 
curve  represents  the  energy  absorbed  by  the  surface/structure. 

Four  idealized  resistance  functions  were  selected  for  analysis 
to  represent  the  primary  resistance  features  illustrated  in 
Figure  4.  In  each  case  no  inertial  effect  of  the  surface/structure 
was  included.  The  four  idealized  resistance  functions  are  illus- 
trated in  Figure  5.  They  include  the  perfect  absorber  which  is 
characterized  by  a constant  resistance  force.  Note  that  very  few 
structures  of  this  nature  will  exist  and  care  must  be  used  to  ensure 
that  adequate  displacement  capacity  (x  ) is  available  so  that  the 
system  will  not  "bottom  out".  The  second  case  is  the  "classical" 
elastic  case.  The  third  and  fourth  cases  treat  the  assumed  to  be 
more  common  cases  of  strain  hardening  yielding  structures.  In  the 
latter  case  an  initial  threshold  force  level  exists. 

An  analysis  of  the  response  of  each  of  the  resistance  func- 
tions can  now  be  performed  to  determine  the  acceleration  history  of 
the  body  component  during  the  impact  process.  This  acceleration 
history  can  then  be  used  to  evaluate  an  appropriate  injury  index 
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Load 


FIGURE  4.  POTENTIAL  LOAD-DEFLECTION  RELATIONSHIP 
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FIGURE  5.  SELECTED  RESISTANCE  FUNCTIONS 
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such  as  the  severity  index,  SI,  or  the  peak  acceleration,  aQ.  For 
the  former,  the  load  application  time  t^  is  also  a factor  both  with 
respect  to  evaluating  the  magnitude  of  the  severity  and  perhaps  in 
defining  the  nature  of  the  injury. 

3.2  Perfect  Absorber 


The  perfect  absorber,  case  A,  is  characterized  by  a constant 

resistance  force,  F . This  force  is  derived  from  the  product  of  the 

application  area,  A , and  the  yield  stress,  a , which  is  a property 

c y 

of  the  material  involved.  Thus 


The  equation  of  motion  for  the  impacting  body  component  of  weight, 
Wc,  is 


dV  F_ 
dt  = " 

c 


(4) 


where 

V = velocity 
t = time 
F = force  = FQ 

Subject  to  the  initial  condition 
V( 0)  = VQ 

and  (5) 

x( 0 ) = 0 


where  x is  the  displacement  associated  with  the  contact  point. 
Integration  yields 


V = 


V. 


V 


(6) 


where  a = F /W  is  the  deceleration  (a  constant).  The  body  compo- 
~ 0 c c 

nent  comes  to  rest  (i.e.,  V=0)  at 


t 


t = V /an 
o 00 


(7) 
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The  displacement  is  given  by  the  differential  equation 


which  when  integrated  yields 

x = V t - a^  t2/2 
oo 


(8) 

(9) 


The  maximum  displacement,  x , is 


x = 1/2  V t 
o oo 


(10) 


The  motion  is  fully  arrested  at  this  point  and  all  of  the  energy 
inherent  in  the  incident  body  component  is  totally  dissipated.  The 
duration  of  the  load,  t^,  is  thus  equal  to 


Vc 

Vy 


(ii) 


Finally,  the  severity  index  (SI)  from  its  basic  definition  (see 
equation  ( 1 ) ) is 


SI  = (a  ) 1 ^ Vn 
o o 


1 

Wcj 

Q 

<< 

1 

1.5 


V 


(12) 


If  the  peak  acceleration  aQ  is  used  as  a measure  of  injury  then  for 
the  perfect  absorber 


* 

= a 

o o 


V* 

w_ 


(13) 


3*3  Elastic  Response 

The  elastic  response  surface,  case  B,  is  characterized  by  a 
resistance  . force  which  is  proportional  to  the  contact  point 
displacement,  that  is 

' F = kx  (14) 

where  k is  the  so-called  spring  constant.  The  above  resistance 
form  is  absolute  and  independent  of  the  direction  of  the  movement 
(i.e.,  loading  versus  unloading  phase).  The  solution  for  this 
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resistance  case,  for  the  above  cited  boundary  conditions,  is 
classic  and  given  for  the  velocity  as 


V(t)  = VQ  cos(wt)  (15) 

where  

-Vh 

The  acceleration  and  displacement  .histories  are  given  as 
a = ajVQ  sin(u)t) 

and  (16) 

x = xQ  sin(ajt) 

where  the  maximum  displacement  is 


The  time  at  which  the  maximum  displacement  is  achieved  is 


Due  to  the'  conservative  nature  of  an  elastic  system  the  energy 
inherent  in  the  motion  of  the  body  component  is  transferred  without 
loss  into  potential  energy  of  deformation  of  the  elastic  system  and 
this  energy  is  readily  available  to  accelerate  the  body  component 
back  in  the  direction  from  which  it  came.  This  response  feature 
results  in  further  acceleration  of  the  body  component  until  such  a 
time  that  it  loses  contact  with  the  surface.  The  above  response 
equations  apply  and  the  load  duration  t^  is  simply 


t 


d 


(19) 


The  severity  index  for  this  case  is 


SI  = 2B. 


0.75 


V 


2.5 
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where 


(20) 


TT  / 2 


0 


The  peak  acceleration  is 


(21) 


3.4  Inelastic  Response 


The  inelastic  response  surface,  cases  C and  D,  are 
characterized  by  a response  force  which  is  similar  in  form  to  that 
of  the  elastic  case  (equation  (14))  except  that  all  of  the  trans- 
ferred energy  is  lost  due  to  plastic  deformation  of  the  surface 
(structure)  and  thus  no  rebound  can  occur.  Equation  (1<4)  only 
applies  for  the  loading  phase. 

Case  C,  the  strain  hardening  case,  corresponds  identically  to 
the  resistance  function  form  of  the  elastic  response  surface  and 
the  motion  solution  equations  (15)  and  (16)  apply.  However,  the 
body  component  comes  to  rest  at  the  point  of  maximum  displacement, 
hence  the  duration  of  the  load  is 


The  peak  acceleration  is  given  by  equation  (21). 

1 Case  D is  the  more  general  case  and  includes  case  A as  one 
limiting  case  and  case  C as  the  other  limiting  case.  The  resistance 
for  case  D includes  an  initial  resistance  Fc  at  zero  displacement, 
thus  the  resistance  is 


(22) 


All  of  the  other  motion  results  apply  directly. 


The  severity  index  for  case  C is  thus 


(23) 


F = F + kX. 
c 


(24) 
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The  corresponding  motion  solutions  are 


coF. 


V = VQ  cos  (tot)  - — — sin(wt) 


a = 


2 Fc 

-wV  sin(ojt)  - to  — cos(tot) 


‘and 


x = 


Vn  F Fn 

— sin(oot)  + t—  cos(wt)  - rr-' 


to 


(25) 


The  maximum  deflection  and  the  load  duration  can  be  readily 
evaluated  by  considering  the  general  case  to  be  composed  of  two 
similar  cases  of  the  form  of  case  C that  is,  the  difference  of  two 
triangular  resistance  functions.  A hypothetical  preload  deforma- 
tion x„  is  used  such  that 
c 


(26) 


The  initial  velocity  for  the  preload  contribution  must  be 
correspondingly  increased  such  that  when  x=0  (and  F=FQ),  V=VQ. 

The  load  duration  is 


t . = t = — arc  tan 
d o to 


and  the  maximum  displacement  is 
V 

xQ  = — sin(oJt0)  + xc  cos(o3t0)  - xQ , 

The  severity  index  (SI)  is 

*|0. 75 


V 

_ / Wc 

nr  v. 

. 3 ' 
o o 
• \x 

1 

= V — arc  tan 

V — — 

v k * F„ 

c 

. (27) 


(28) 


SI  = f(v) 


k_ 

W. 


2.5 


(29) 


where 

F(v) 


B(  3) 

(cos  8) 2*  ^ 


B(6) 


7T/2 

f . 2. 

= sm 


ad  a 
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3 = arc  tanv 


v = ps=‘ 

The  function  f(v)  has  been  evaluated  and  is  presented  in  Figure  6. 
-In  the  limit  for  large  values  of  v 

f(v)  = v1,5.  (30) 

The  peak  acceleration  is  given  as 


3*5  Sensitivity  Evaluation 

An  examination  of  the  dependence  of  the  severity  index  and  the 
peak  acceleration  on  the  parameters  involved  will  be  of  interest 
inasmuch  as  such  a sensitivity  evaluation  will  indicate  which  of 
the  parameters  are  most  significant  in  this  injury  analysis.  In 
this  manner  appropriate  attention  can  be  given  to  those  parameters 
in  the  analyses  presented  in  the  subsequent  chapters  of  this  report 
and  in  the  allocation  of  resources  for  further  research. 

The  parameters  which  are  involved,  although  not  necessarily  in 
all  cases,  are  the  impact  velocity,  VQ,  the  critical  force,  Fq,  the 
spring  constant,  k,  and  the  weight  of  the  body  component,  WQ.  The 
severity  index  (SI)  equations  can  be  expressed  in  a logarithmic 
differential  form  as 


d(SI) 

SI 


dV  dF 

C — + r 9. 

V + LF 
0 c 


(32) 


where  the  coefficients  C.  are  defined  as  influence  coefficients 

1 «' 

and  correspond  to  the  contribution  that  a percentage  change  in  the 
ith  parameter  makes  to  the  percentage  change  in  the  severity  index. 
Appropriate  differentiation  of  the  severity  index  equations 
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u 


FIGURE  6.  SEVERITY  INDEX  FUNCTION  f (v)‘ 
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presented  above  provide  the  values  of  these  coefficients  for  each 
of  the  four  response  cases  treated.  A similar  treatment  was 
applied  to  the  peak  acceleration  solutions  to  yield  corresponding 
influence  coefficients  K^,  that  is 


d(a*) 


= K 


dVo  dFc  dk  dWc 

V V F,  ^k  V W 


(33) 


The  values  of  the  influence  coefficients  are  summarized  in  Table  1. 
For  the  general  case,  case  D,  the  coefficients  are  functions  of  the 
parameter  v and  must  be  treated  as  infinitesimal  influence  coeffi- 
cients. For  this  general  case  the  values  of  the  influence 
coefficients  for  the  severity  index  and  peak  acceleration  are 
presented  in  Figures  7 and  8 respectively. 


TABLE  1. -SUMMARY  OF  INFLUENCE  COEFFICIENTS 

Case  (SI)  a* 

o 


CV 

CF 

Ck 

°W 

KV 

kf 

Kk 

A 

1.0 

1.5 

0 

-1.5 

0 

1.0 

0 

-1.0 

B 

2.5 

0 

0.75 

-0.75 

1.0 

0 

0.5 

-0.5 

C 

2.5 

0 

0.75 

-0.75 

1.0 

0 

0.5 

•0.5 

D* 

functions 

of  V 

functions 

of  V 

see  Figure  7 and  8 


An  examination  of  the  absolute  values  of  the  influence  coeffi- 
cients indicates  that  the  single  most  important  parameter  is  the 
impact  velocity,  VQ.  Clearly,  then  the  analysis  of  the  motion  of 
passengers  in  the  various  passenger  configurations  must  be  designed 
to  evaluate  the  impact  velocity  accurately  and  to  identify  the 
factors  which  influence  its  value.  The  weight  of  the  body 
component,  W , is  moderately  important  but  it  does  not  represent  a 
parameter  which  could  be  controlled  by  car  interior  component 
considerations . 
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FIGURE  7.  INFLUENCE  COEFFICIENTS  FOR  SEVERITY  INDEX 


Case 
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FIGURE  8.  INFLUENCE  COEFFICIENTS  FOR 


Similarly,  the  application  area,  A , is  largely  a "people"  param- 
eter.  Since  the  force,  for  case  A,  is  proportional  to  the  yield 
stress  the  influence  coefficient  for  the  stress,  a , is  essentially 
equal  to  that  of  the  force.  The  value  of  the  influence; coefficient 
for  force  is  moderate.  The  least  significant  parameter  appears  to 
be  the  spring  constant  or  stiffness  parameter.  Although  this 
“parameter  may  be  the  least  influential  it  is  one  with  which  the 
greatest  improvements  may  be  achieved  by  design  because  of  the 
greater  percentage  change  that  may  be  available. 

3.6  Composite  Systems 

The  examined  resistance  functions  were  quite  simple  and  do  not 
represent  real  systems  with  great  precision.  The  resistance  of 
real  structural  systems  to  local  impact  will  be  much  more  complex 
and  depend  upon  the  point  of  impact,  inertial  effects  of  the 
surface  or  structure,  and  many  other  factors.  In  this  section  the 
variation  of  the  stiffness  of  a simple  composite  beam/support 
system  is  examined.  Such  a system  could  approximate  some  typical 
car  interior  components  such  as  partitions,  stanchions,  and  perhaps 
to  a lesser  degree  some  seatbacks.  The  beam/support  system  consid- 
ered here  is  a uniform  simply  supported  beam  mounted  upon  two 
identical  spring  supports. 

Consider  first  the  two  limiting  cases  of  a flexible  beam  on 
rigid  supports  and  that  of  a rigid  beam  on  flexible  supports.  The 
spring  constants  A for  a load  applied  at  a position  A along  the 
beam  (A  = 0 corresponds  to  the  support  position  and  A = 0.5 
corresponds  to  the  midspan  position)  is  presented  in  a normalized 
form  in  Figure  9(a).  The  spring  constant  has  been  normalized  by 
the  minimum  values  respectively  for  each  case.  The  stiffness 
variation  for  different  application  points  is  reasonably  constant 
for  the  rigid  beam  on  spring  supports  with  a maximum  variation  of  a 
factor  of  two.  On  the  other  hand,  the  system  of  the  flexible  beam 
on  rigid  support  results  in  a very  large  variation  for  different 
application  points. 
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FIGURE  9.  LOCAL  STIFFNESS  FOR  COMPOSITE  BEAM/SUPPORT  SYSTEM 
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In  particular,  the  system  becomes  very  stiff  near  the  support  posi- 
tions. Approximately  20  percent  of  the  possible  impact  positions 
result  in  local  stiffnesses  which  are  at  least  10  times  greater 
than  the  miminum  stiffness. 

The  design  of  the  overall  system  can  be  improved  with  respect 
-to  body  component  impact  considerations  by  distributing  the  stiff- 
ness between  the  two  components  of  the  system.  The  local  stiffness 
for  this  general  case  is  presented  in  Figure  9(b).  This  system 
stiffness  is  normalized  by  the  minimum  or  midpoint  stiffness  of 
the  beam  component.  A much  more  balanced  design  can  thus  be 
achieved.  In  particular,  the  extremely  high  local  stiffness  near 
the  support  position  can  be  significantly  reduced. 
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Yield  stress,  a (psi) 

J « 

FIGURE  10.  INFLUENCE  OF  YIELD  STRESS  ON  SEVERITY  INDEX  - 
PERFECT  ABSORBER 
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3.7  Typical  Results 


The  severity  index  was  evaluated  for  the  case  of  a perfect 
absorber  and  the  case  of  an  elastic  response.  Recall  that  the 
latter  case  also  corresponds  to  case  C by  reducing  the  severity 
index  by  a factor  of  two.  The  head  impact  situation  was  selected 
-for  this  evaluation  and  nominal  values  for  the  body  component 

weight  and  application  area  of  10  lb  and  10  sq  inches  were 

selected.  The  results  cover  the  impact  velocity  range  of  from  10  to 
50  f ps . The  results  for  the  perfect  absorber  are  presented  in 

Figure  10  and  suggest  that  a yield  stress  in  the  range  of  from  80  to 
200  psi  would  produce  a moderate  injury  for  the  velocity  range 

cited.  The  results  for  the  elastic  response  case  are  presented  in 
Figure  11  and  illustrate  the  large  role  that  the  impact  velocity 
plays  in  this  system.  A value  of  a spring  constant  of  1000  lb/inch 
can  produce  rather  severe  injuries  at  all  but  the  very  low  velocity 
levels . 

3.8  Nomenclature 

Symbols  and  definitions  used  in  Section  3: 
a acceleration 


Bo 

ci 

CF 

C* 

CV 

CW 

f 

F 


specific  acceleration 
peak  acceleration 
application  area 
a function 
a constant 

ith  influence  coefficient 

severity  index  influence 

severity  index  influence 

severity  index  influence 

severity  index  influence 

a function 
force 


for  severity  index 


coefficient 

for 

force 

coefficient 

for 

spring  constant 

coefficient 

for 

velocity,. 

coefficient 

for 

weight 
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FIGURE  11.  INFLUENCE  OF  SPRING  CONSTANT  ON  SEVERITY  INDEX 
ELASTIC  BEHAVIOR 


ri 

X.v. 


K, 


K 


V 


SI 

t 

*0 

V 

Vo 

x 

x_ 


Wc 

a 

6 

X; 

v 

Cl) 

A 

A 


critical  force 
maximum  force 
spring  constant 

ith  influence  coefficient  for  peak  acceleration 

peak  acceleration  influence  coefficient  for  force 

peak  acceleration  influence  coefficient  for  spring  constant 

peak  acceleration  influence  coefficient  for  velocity  K Deak 

acceleration  influence  coefficient  for  weight 

severity  index 

time 

load  duration 

time  of  maximum  displacement 
velocity 

incident  velocity 
displacement 

hypothetical  displacement 

maximum  displacement 

generic  displacement 

generic  displacement 

generic  displacement 

body  component  weight 

dummy  variable  for  integration 

dimensionless  variable 

dimensionless  position  along  beam 

dimensionless  variable 

frequency 

normalized  stiffness 
relative  normalized  stiffness 
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4.  PASSENGER  MOTION 


The  primary  goal  in  this  section  is  to  develop  motion  details 
for  persons  in  the  four  basic  passenger  configurations,  and  to 
identify  the  nature  and  magnitude  of  impact  velocities  of  their 
equivalence  in  each  of  these  situations.  As  a prelude  to  present- 
ing these  results  it  was  desirable  to  examine  a class  of  train  car 
acceleration  pulses  to  establish  their  characteristics  and  ulti- 
mately to  unify  these  pulses  by  developing  equivalence  parameters. 
A simple  lumped  mass  translation  behavior  analysis  was  then 
performed  to  obtain  an  insight  into  potential  motion  environment. 
Finally  a brief  summary  of  anthropometric  data  is  presented  to 
establish  some  of  the  needed  passenger  characteristics. 


At  this  point  the  specific  motion  details  of  passengers  in  the 
standing,  forward  facing  seated,  backward  facing  seated,  and  side- 
ward facing  seated  passenger  configurations  are  presented.  The 
section  ends  with  a brief  and  somewhat  superficial  treatment  of 
multiple  passenger  interactions.  This  last  effort  represents  an 
attempt  to  "measure”  the  potential  influence  of  some  real  world 
situations . 

4.1  Acceleration  Pulse  Shape  Characteristics 

An  examination  of  some  numerical  results  obtained  during  the 
course  of  the  current  program  dealing  with  the  motion  details  of 
railcars  during  a crash,  indicate  that  the  acceleration  histories 
can  be  idealized  with  relatively  simple  pulse  shapes.  These 
include  the  step  pulse,  a double  step  pulse,  and  a triangular  pulse 
shape.  Some  typical  acceleration  histories  are  illustrated  in 
Figure  12(a)  and  12(b).  It  is  clear  that  some  fine  structure  will 
exist  and  in  addition  that  some  low  amplitude  details  may  persist 
for  some  time.  It  appears  reasonable  to  define  several  idealized 
acceleration  pulse  shapes  and  to  attempt  to  utilize  these  in  the 
subsequent  passenger  motion  studies. 

Two  rather  general  pulse  shapes  are  selected,  one  is  a trape- 
zoidal pulse,  the  other  a double  step  pulse.  These  pulse  shapes  are 
illustrated  in  Figures  12(c)  and  12(d). 
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FIGURE  12.  ACCELERATION  PULSE  CLASSIFICATION 
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The  trapezoidal  pulse  includes  within  its  framework  the  basic  step 
pulse  and  the  full  spectrum  of  triangular  pulses.  The  double  step 
pulse  also  includes  with  its  form  the  basic  step  pulse.  It  appears 
reasonable  and  practical  to  refer  all  pulse  shapes  to  the  basic 
step  pulse;  a pulse  shape  which  itself  is  a good  approximation  to 
many  observed  (i.e.,  calculated)  acceleration  pulses. 

By  definition  the  area  under  the  acceleration  pulse  a(t)  is 
equal  to  the  change  in  velocity  VQ  of  the  train  car 


V 


o 


o 

a(t)  dt 


(34) 


where  t = time  and  t = pulse  duration.  It  follows  from  these 

o 

definitions  that  the  average  acceleration,  a , is  given  as 


Two  parameters  have  been  elected  as  useful  parameters  with  which 
to  characterize  these  acceleration  pulses.  One  is  an  acceleration 

peaked  factor,  y,  which  is  simply  the  ratio  of  the  peak  accelera- 

* 

tion  aQ  to  the  average  acceleration,  that  is 


Y 


(36) 


The  other  is  a stopping  distance  factor  3 defined  as 

3 = 


Vo 


(37) 


where  xQ  = stopping  distance. 


The  motion  details  corresponding  to  a basic  step  pulse  accel- 
eration, i.e.,  a constant  deceleration  level  a , are  defined  by 

a(t)  = aQ  0 < t < tQ 

V = V -at  (38) 

o o 


X = V - 


V 
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where  V 

= velocity 

X 

= distance  or  displacement. 

When  the 

train  car  comes  to  rest  (V  = 0)  and 

t = 

fco  ; 

x = 

X 

0 

“hence  one  obtains  equation  (35)  from  the  velocity  time  relationship 


and 

xo 

1/2  VQt0  (39) 

from  the 

displacement  relationship.  Thus  for  the  basic  step  pulse 

Y = 

1 

3 = 

1/2. 

These  two  parameters  were  evaluated  for  both  the  trapezoidal  and 
the  double  step  pulse  shape.  For  the  trapezoidal  pulse  one  obtains 


Y = 

(40) 

1 + VT1 

and 

Q 

( 1-t2)  (3-t2)2  + t2(3-t22)  - 2 

p = 

3t1-T2;(i  + T2-V 

where,  from  Figure  12(c) 


T1  ; 

. ll  < T 

' - 2 

and 

-r 

t2 

T2  : 

= V 

These  results  are  presented  in  Figure  13*  Note  that  the  stopping 
distance  factor,  3,  only  varies  from  one-third  to  two-thirds  with 
most  of  the  pulse  domain  having  a value  of  3 of  approximately  one- 
half.  The  acceleration  peaked  factor,  Y , varies  from-- one  for  the 
basic  step  pulse  to  two  for  the  entire  spectrum  of  triangular 
pulses . 
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FIGURE  13.  CHARACTERISTICS  OF  TRAPEZOIDAL  ACCELERATION  PULSE 
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Similar  results  were  obtained  for  the  double  step  pulse  shape, 
namely 


y ~ 


1 


1 


5 - 


T16 


(41) 


and 


0. 5(6-5  xi  + 
6 = Tc 


1 + 6 - T-j  6) 


where  from  Figure  12(d) 
t 


1 


r 


and 


<5  = -T  * 


These  results  are  presented  in  Figure  14.  If  one  excludes  the  lower 
left-hand  position  of  the  pulse  domain  on  the  basis  that  any  re- 
shape pulse  could  be  better  approximated  by  a shorter  duration 
basic  step  pulse  then  it  is  clear  that  the  range  over  which  both 
factors  vary  is  not  excessive.  This  observation  suggests  that  the 
basic  step  pulse  can  be  used  for  analysis  purposes  in  many  cases  and 
that  adjustment  can  be  made  for  other  pulse  shapes  by  using  the 
above  factors  in  some  prescribed  manner  depending  upon  the  kind  of 
equivalent  required.  As  an  example,  for  energy  (or  work)  equiva- 


lence the  kinetic  energy  0.5  would  be  conserved,  that  is 

9 V ax 

0.5  = 0.5  — (VQt0)  = — 2g—  = constant  - 

o p 


(42) 


Thus  the  proper  relationship  between  the  reference  and  equiva- 
lent system  requires 


(aoxo) 


equiv  = 


'Vo5 


ref 


(43) 


The  four  motion  parameters,  V , aQ,  xq  and  tQ  are  related  by 
two  equations  (equations  (35)  and  (37)  for  any  given  pulse  shapq). 
Thus  one  can  conveniently  construct  a solution  diagram  for  any 
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FIGURE  14.  PULSE  CHARACTERISTICS  FOR  THE  DOUBLE  STEP  PULSE 
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pulse  shape  in  which  the  specification  on  any  two  parameters  yields 
the  other  two.  Two  such  diagrams  were  constructed  and  are  pre- 
sented in  Figures  15  and  16.  The  first  applies  to  the  step  pulse 
and  provides  a perspective  of  the  impact  domain  of  general  interest 
to  the  intracity  train  crash  problem.  The  second  applies  to  one 
particular  triangular  pulse  as  well  as  any  pulse  shape  with  a 
“stopping  distance  factor  equal  to  two-thirds. 

4.2  ’’Lumped  Mass”  Translation  Behavior 

Perhaps  the  simplest  representation  of  the  motion  of  a passen- 
ger within  a compartment  of  a train  which  is  abruptly  decelerated, 

as  in  a crash,  is  that  of  a lumped  mass  initially  located  some 

* 

distance,  x , from  an  obstacle  such  as  a seatback  or  bulkhead, 
fixed  to  the  compartment,  and  moving  with  the  train.  When  the  crash 
occurs  at  zero  time  (t=0)  the  compartment  (i.e.,  bulkhead)  moves  in 
accordance  with  some  prescribed  acceleration  pulse  as  previously 
mentioned.  The  passenger  originally  moving  at  the  train  velocity, 
V , continues  to  move  until  impact  with  the  bulkhead  occurs.  The 
train  velocity,  V , corresponds  to  the  change  in  velocity  experi- 
enced by  the  car  during  a crash  and  may  not  be  equal  to  the  actual 
velocity  of  the  car.  The  impact  velocity  that  is  the  relative 
velocity  between  the  passenger  and  the  bulkhead  is  of  interest 
because  it  is  a major  parameter  in  the  production  of  injury 
mechanisms.  It  appears  reasonable  with  this  simple  model  to 
introduce  a nominal  average  force  (expressed  as  a deceleration)  as 
an  additional  parameter  to  examine  the  influence  of  a variety  of 
possible  motion  retarding  mechanisms,  such  as  seat  friction,  may 
have  upon  the  impact  velocity. 

Representative  paths  are  illustrated  in  the  path  diagram  of 

Figure  17.  The  following  analysis  is  presented  in  a dimensionless 

f^erm  to  facilitate  its  interpretation  and  to  expand  its  usefulness. 

The  time  scale  has  been  made  dimensionless  by  normalising  time  by 

the  crash  duration,  t . The  position  scale  has  been  made 

o 

dimensionless  by  dividing  the  absolute  position  x by  the  stopping 

distance,  x . thus 
7 o 
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FIGURE  15.  SOLUTION  DIAGRAM  FOR  STEP  PULSE 
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FIGURE  16.  SOLUTION  DIAGRAM  FOR  TRIANGULAR  PULSE 
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FIGURE  17.  PATH  DIAGRAM  FOR  PASSENGER/BULKHEAD  INTERACTION 
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Zero  time  is  defined  as  the  start  of  the  crash  period.  Zero  posi- 
tion is  defined  as  the  location  of  the  train  bulkhead  at  zero  time 
and  positive  positions  increase  in  a direction  away  from  the 
bulkhead.  The  bulkhead  path  illustrated  in  Figure  17  corresponds 
to  a step  acceleration  pulse.  The  bulkhead  motion  vanishes,  by 
definition  at  £ = 1 and  x = 1.  Thereafter  it  is  motionless.  At 
zero  time  the  bulkhead  and  all  of  the  passengers  (P^  and  P^)  move  at 

a velocity  V in  the  negative  c direction.  The  passenger  is 

0 * 

initially  located  at  l,  = c where 


The  passenger  path  in  Figure  17  labeled  B corresponds  to  a passen- 
ger originally  at  P^  who  moves  unretarded  with  the  velocity  VQ  and 
impact  occurs  with  the  bulkhead  at  point  IB.  The  bulkhead  is  still 
in  motion  hence,  the  impact  velocity  is  somewhat  less  than  VQ. 
The  passenger  path  labeled  A corresponds  to  a passenger  originally 
at  who  also  moves  in  an  unretarded  fashion.  Impact  with  the 
bulkhead  occurs  at  point  IA  at  which  time  the  bulkhead  velocity  is 
zero.  Therefore,  the  impact  velocity  is  VQ.  Finally,  the  passen- 
ger path  labeled  C corresponds  to  the  passenger  originally  located 
at  P^,  however  in  this  case  the  passenger  motion  is  retarded  such 
that  impact  with  the  bulkhead  occurs  at  point  IC.  The  bulkhead 
velocity  is  zero  at  this  time,  however,  the  velocity  of  passenger 
is  greatly  reduced  such  that  the  impact  velocity  is  considerably 

less  than  V . Recall  that  the  injury  level  may  be  proportional  to 
"■pc  O 

V:  . The  impact  velocity,  Vi , is  normalized  by  the  initial 

velocity,  V , to  provide  a dimensionless  impact  variable,  viz. 


x 

X 


(45) 


o 


(46) 
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and  the  retarding  acceleration  a’  is  normalized  by  the  average 
acceleration  of  the  train  car,  a , to  form  a dimensionless  retard- 
ing force,  viz. 


(47) 


-Two  impact  domains  exist;  one  corresponds  to  impacts  which  occur 
while  the  bulkhead  is  still  in  motion  and  the  other  corresponds  to 
impacts  which  occur  after  the  bulkhead  motion  stops.  This  latter 
domain  has  been  defined  as  the  far  field  domain.  The  impact  veloc- 
ity in  this  domain  is 

ii  = l/1  - 2a6(  5*+1)  W) 

and  applies  to  any  pulse  shape.  The  parameter  3 is  the  stopping 
distance  factor  defined  in  the  previous  section.  The  far  field 

results  for  the  step  pulse  case  are  presented  in  Figure  18.  Recall 

* / \ 

that  £ is  the  (normalized)  location  of  the  passenger  at  the  time 
of  the  crash.  For  this  pulse  shape  8 = 0.5,  therefore  equation 
(48)  becomes 


i = “\Jl  ~ ot(5*+1)  • 


(49) 


The  results  of  Figure  18  are  generally  applicable  to  other  pulse 
shapes  if  the  value  of  a used  in  this  figure  is  viewed  as  an  equiva- 
lent a equal  to  the  product  2a3  for  other  pulse  shapes.  The  region 
of  applicability  of  this  far  field  relationship  depends  upon  the 
shape  of  the  acceleration  pulse.  Thus  it  will  be  desirable  to 
examine  some  near  field  results. 


The  near  field  results  for  the  step  pulse  are  shown  in 
Figure  19*  The  near  field  region  for  this  pulse  shape  is  bounded  by 
the  curves  a = 0 and  a = a The  latter  is  given  as 

a = 1 - ;*  (50) 

along  which 

«*i  = C • 
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FIGURE  18.  IMPACT  ENVIRONMENT  - LARGE 


n 

a(t) 


FIGURE  19.  IMPACT  ENVIRONMENT  - STEP  PULSE 


48 


Along  the  former 


These  results  illustrate  that  in  the  far  field,  reductions  in 
impact  velocity  can  be  achieved  by  a combination  of  increased 
-s-'eparation  distance  and  increased  retardant  force.  In  the  near 
field  the  impact  velocity  can  be  reduced  by  reducing  the  separation 
distance  and/or  increasing  the  retardant  force.  It  appears  sig- 
nificant to  identify  the  impact  zone  which  applies,  or  its  approxi- 
mate equivalence,  the  stopping  distance,  xQ . 

The  results  become  more  difficult  to  interpret  when  one 
examines  other  pulse  shapes.  One  of  the  right  triangular  pulses 
was  used  in  a similar  analysis.  These  results  are  presented  in 
Figure  20.  The  far  field  results  are  unchanged,  as  was  indicated 
above,  however,  the  domain  of  application  has  changed  somewhat. 
The  role  of  the  retarding  force  becomes  more  obscure.  In  effect, 
what  can  happen  is  that  an  increase  in  the  retarding  force  results 
in  a delay  in  the  impact  time  and  due  to  larger  late  time  decelera- 
tion associated  with  this  particular  pulse  shape  the  slightly 
delayed  impact  is  more  severe  than  an  earlier  impact  would  have 

been.  This  is  seen  more  clearly  in  Figure  21  where  the  influence  of 

* 

a is  presented  along  lines  of  constant  z,  , for  both  pulse  shapes. 

Other  pulse  shapes  (within  the  trapezoidal  category)  were 
examined  and  similar  apparent  anomalies  are  present.  Their 
exploration  at  this  time  will  not  be  particularly  productive.  It 
may  be  useful  to  define  the  boundary  between  the  near  field  and  far 
field  for  the  general  trapezoidal  pulse  shape,  that  is  the  ^ = a 

v 

path.  This  is  given,  in  general  as 

«c  = 2 - 2g(1  + £*) 

and  (52) 

ii  --  1 - ac 
along  this  path. 
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FIGURE  20.  IMPACT  ENVIRONMENT  - TRIANGULAR  PULSE 
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FIGURE  21.  INFLUENCE  OF  a ON  IMPACT  VELOCITY 


4.3  Passenger  Characteristics  - Anthropometric  Data 

The  tumbling  motion  that  a human  being  experiences  from  a 
given  motion  stimulus  will  depend,  in  part,  upon  various  anthropo- 
metric parameters.  Furthermore,  the  values  of  these  parameters 
will  generally  be  statistical  in  nature,  although  some  parameters 
-may  vary  in  a systematic  manner  with  age,  sex,  and  other  classifi- 
cation factors.  For  the  purpose  of  this  study  a number  of  parameter 

values  and  their  variations  must  be  established  or  estimated. 

4 

Damon,  et  aL,  was  the  source  of  most  anthropometric  data  used  in  the 
subsequent  analysis.  In  some  instances  simple  approximate  scaling 
relations  were  established  to  encompass  parameter  values  in  a 
systematic  manner  for  certain  population  classifications. 

The  general  type  of  anthropometric  data  that  is  required  to 
evaluate  tumbling  motions  of  passengers  are  whole  body  or  body 
component  weight,  one  or  more  characteristic  lengths  such  as  height 
and  location  of  the  body  c.g.  relative  to  some  reference  point,  and 
the  body  or  body  component  moment  of  inertia  about  some  specific 
axis  passing  through  the  appropriate  c.g.  These  data  will 
generally  be  presented  as  mean  or  average  values  and  may  be 
supplemented  with  some  indication  of  their  variance  such  as  their 
standard  deviations.  Unfortunately,  anthropometric  data  are  quite 
limited  and  are  frequently  biased  by  the  use  of  rather  similar  and 
limited  population  subgroups,  such  as  young-white-adult-males. 
Nevertheless,  reasonable  average  values  can  be  estimated  together 
with  somewhat  less  reliable  variance  estimates.  These  anthropo- 
metric uncertainties  must  be  combined  with  other  types  of 
uncertainties  and  will  result  in  a statistical  smearing  of  the 
final  results. 


Design”,  Harvard 


Damon,  A.,  et  al.,  "The  Human  Body  in  Equipment 
University  Press,  1966. 
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To  demonstrate  the  general  nature  of  the  anthropometric  data 
and  to  establish  scaling  relationships  for  use  in  the  analysis  of 
the  motion  characteristics  of  the  standing  passenger  height/weight 
data  are  presented  in  Figure  22.  The  average  adult  ;male  weighs 
1'65  lb  and  is  5.77  ft  (5  ft  9.25  inches)  tall.  Height  distribution 
data  suggest  limiting  values  which  can  be  interpreted  as  "tail"  and 
"short"  category  classifications.  These  values  are  illustrated  in 
Figure  22'  and  are  listed  in  Table  2,  and  represent  approximately 
95  percent  of  the  male  adult  group.  Similarly,  but  from  a smaller 
data  base,  the  weight  variance  data  permit  an  estimate  of  corre- 
sponding limiting  values  for  the  "thin"  and  "fat"  category  classi- 
fication. These  values  are  also  illustrated  in  Figure  22  and 
listed  in  Table  2.  The  weight  variance  is  referred  to  herein  as  a 
"shape"  parameter.  Thus  eight  shape-height  categories  can  be 
established  in  addition  to  the  average  adult  male  (now  identified 
by  the  avg-avg  notation).  These  categories  bound  the  adult  male 
domain  as  illustrated  in  Figure  22  and  include  more  than  95  percent 
of  the  corresponding  population. 


TABLE  2. -HEIGHT/WEIGHT  CHARACTERISTICS 


Category 

(shape-height) 

Weight 

(lb) 

Height 

(ft) 

Moment 
Inertia 
(lb  sec^ 

of 

•kk 

ft) 

X 

a 

Avg-Avg  * 

165 

5. 

77 

8.58 

1. 

.00 

1, 

.00 

Avg- Short 

124 

5. 

25 

5.88 

0. 

.91 

1, 

.00 

Avg-Tall 

210 

6. 

25 

11.81 

1. 

.08 

1. 

.00 

Thin- Short 

90 

5. 

25 

5.00 

0. 

,91 

0. 

.85 

Thin- Avg 

119 

5. 

77 

7.29 

1. 

.00 

0. 

.85 

Thin-Tall 

152 

6. 

25 

10.04 

1. 

,08 

0. 

,85 

Fat-Short 

165 

5. 

25 

6.76 

0. 

.91 

1. 

.15 

Fat -Avg 

218 

5. 

77 

9.87 

1. 

.00 

1. 

,15 

Jat-Tall 

278 

6. 

25 

13.58 

1. 

.08 

1. 

,15 

"10-yr-old  child 

78 

4. 

50 

3.17 

0. 

.78 

1. 

,00 

6-yr-old  child 

40 

3. 

60 

1.30 

0. 

62 

1. 

00 

Adult  female  is 

•kit 

Standing 

similar  to 

avg- 

short . 
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FIGURE  22.  HEIGHT/WEIGHT  DATA  SUMMARY 


Height  weight  data  for  both  males  and  females  of  various  age 
groups  from  birth  to  age  19  years  is  shown  in  summary  fashion  in 
Figure  22.  These  data  are  consistent  with  the  average  adult  male 
data  and  suggest  that  the  average  adult  female  category  may  be 
defined  by  values  corresponding  approximately  to  the  avg-short 
category  for  adult  males.  The  5-  and  10-year-old  child  categories 
were  also  selected  for  subsequent  analysis  and  the  corresponding 
height/weight  values  are  illustrated  in  Figure  22  and  represented 
in  Table  2. 


Additional  body  component  size,  c.g.  locations,  and  moment  of 
inertia  data  will  be  needed  for  the  subsequent  motion  analyses  and 
are  generally  not  available  for  all  of  the  categories  listed  in 
Table  2.  Thus  several  simple  scaling  laws  were  developed  by 
assuming  that  the  average  body  density  and  internal  density 
distributions  were  similar  to  the  body  shape.  Two  parameters  were 
then  selected  to  represent  the  category  subgroups.  These  are  the 
height  scale  factor  X and  the  shape  scale  factor  a . They  are 
defined  as 


and 


o = 


w 

w 


(53) 


where 

H = height 

H = height  of  average  male  adult 
w = width 

w = width  of  average  male  adult* 

The  corresponding  values  are  presented  in  Table  2.  The  values  of  o 
were  estimated  from  height-weight  data  and  the  appropriate  scaling 
law.  The  scaling  laws  are 
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(54) 

X 

aX 

where 

W = body  weight 

W = body  weight  for  average  male  adult 

I = moment  of  inertia 

T = moment  of  inertia  for  average  male  adult 

L = c.g.  location 

L = c.g.  location  for  average  male  adult 

S = width 

S = width  for  average  male  adult 

"3 

The  height-weight  scaling  law  is  plotted  in  Figure  22  and  illus- 
trates its  appropriation  for  the  "average”  values  from  age  5 to 
adult  age.  These  scaling  laws  were  used  to  generate,  as  needed, 
anthropometric  values.  Moment  of  inertia  estimates  for  the  standing 
passenger  are  listed  in  Table  2.  The  size  and  weight  data  suggest 
that  the  corresponding  variances  are  relatively  constant  for  both 
males  and  females  and  independent  of  age  (above  a few  years  of  age). 

Anthropometric  data  are  also  needed  for  different  body 
positions,  body  movements,  and  body  components.  These  more 
s-pecific  data  are  rather  limited.  Estimates  were  made  and  are 
listed  in  Table  3 for  the  seated  passenger  configurations.  Vari- 
ances (standard  deviation)  of  the  c.g.  location  were  from  5 to 
10  percent,  whereas  variances  of  the  moment  of  inertia  were  15  to 
20  percent.  The  c.g.  location  corresponds  to  the  distance  between 
the  body  component  c.g.  and  the  axis  of  rotation  assumed  in  the 


W 
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L 
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corresponding  analysis.  The  value  of  the  moment  of  inertia 
corresponds  to  rotation  about  the  axis  through  the  body  components 
or  components  of  interest.  The  characteristic  dimension  corre- 
sponds to  the  distance  between  the  potential  impact  location  and 
the  axis  of  rotation  used  in  the  corresponding  analysis.  No 
specific  value  is  needed  for  the  backward  facing  passenger  since 
t-he  injury  mode  is  that  of  whiplash. 


TABLE  3- -ANTHROPOMETRIC  DATA  FOR  SEATED  PASSENGERS 


Seated 

Forward 

Passenger 

Backward 

Facing 

Sideward 

Body  Component 

Torso 

Head 

Whole  Body 

Body  Component  Weight  '(lb) 

87.0 

10.00 

165.00 

Moment  of  Inertia  (lb-ft  ) 

27.4 

0.27 

166.00 

Center  of  Gravity  Location  (ft) 

1.0 

0.30 

0.27 

Characteristic  Dimension  (ft) 

2.5 

2.50 

4.4  Standing  Passenger 

The  response  of  a standing  passenger  to  the  sudden  change  in 
motion  of  the  floor  will  depend  upon  many  factors,  such  as  the 
nature  and  intensity  of  the  change  in  motion  of  the  floor,  the  size 
and  orientation  of  the  passenger,  and  the  passenger’s  alertness  and 
ability  to  respond  to  the  excitation.  If  the  passenger  tumbles  or 
falls  the  passenger  may  strike  the  floor,  nearby  walls  or 
partitions,  or  other  nearby  interior  car  components.  Since  any 
realistic  analysis  will,  of  necessity,  exclude  many  of  the  above 
factors,  the  results  will  have  to  be  interpreted  as  a nominal,  or 
perhaps  upper  limit  experience  of  the  subject  event.  It  is  clear 
that  the  passenger's  use  of  his  or  her  arms  or  hands  may  mitigate 
the  consequences  of  the  event  and  thus  statistically  the  conse- 
quences of  the  event  should  be  reduced. 


IT 

A tumbling  man  (block)  model  was  selected  for  the  analysis  of 
the  standing  passenger.  Although  this  is  a simple  representation 
of  the  passenger  the  use  of  additional  linked  elements  does  not  a 
priori  improve  the  accuracy  of  the  results.  Additional,  usually 
unknown  or  poorly  defined  parameters  are  required  to  properly 
-couple  the  linked  elements  such  as  the  head,  upper  torso,  lower 
torso,  arms  and  legs  together.  The  analysis  of  three-  and  seven- 
element  models  exposed  to  explosion  effects  yielded  essentially 
the  same  results  as  did  the  simple  block  model. 

A series  of  tumbling  block  analyses  were  conducted  for  the 
average  female  subject  to  several  constant  deceleration  floor 
environments  and  velocity  changes.  One  result  is  illustrated  in 
Figure  23(a).  In  the  absence  of  any  nearby  barrier  or  obstacle  the 
standing  passenger  rotates  (head  first)  in  the  direction  of  motion 
and  a head  impact  occurs.  In  the  example  shown  the  train  car  stops 
in  0.4  sec  over  a distance  of  approximately  9 ft.  The  head  impact 
occurs  at  0.66  sec  and  the  impact  velocity  is  approximately 
21.8  fps.  If  a barrier,  such  as  a wall,  was  located  nearby  then  the 
passenger  head  would  have  impacted  at  the  times  and  velocities 
indicated  in  Figure  23(b)  for  the  barrier  distances  indicated.  If 
the  barrier  were  located  at  21.5  ft  for  this  example,  the  top  of  the 
passenger's  head  would  have  impacted  the  barrier  at  a velocity  of 
40.0  fps  simultaneously  with  a floor  impact  of  21.8  fps.  The 
latter  impact  involves  the  side  (or  front,  etc.)  of  the  head.  Such 
multiple  impact  situations  are  unlikely  and  will  not  be  considered 
in  the  analyses. 

An  examination  of  the  motion  of  the  standing  passenger 
revealed  the  following  common  behavior.  Foot  contact  with  the 
floor  was  lost  rather  early  in  the  process  such  that  the  horizontal 
velocity  of  the  c.g.  (&  ) decreased  slightly  relative  to  the 

' C.g. 


5Longinow,  A.,  et  al , "Debris  Motion  and  Injury  Relationships  in 
All  Hazard  Environments",  IIT  Research  Institute,  Final  Report 
J6334,  July  1976. 
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Barrier  Location , x^=5  ft 
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FIGURE  23.  TUMBLING  PASSENGER  TRAJECTORIES 


initial  train  velocity  (i.e.,  the  original  velocity  of  the  passen- 
ger) and  then  remained  constant  at  a value  of  approximately  2 fps 
less  than  the  initial  train  velocity.  The  horizontal  velocity  of 
the  right  side  of  the  head  (x^)  relative  to  the  initial  train 
velocity  increased  initially,  reaching  a maximum  of  3*6  fps  at 
0.3  sec  and  decreasing  to  a velocity  approximately  4 fps  less  than 
the  initial  velocity  at  0.66  sec  when  floor  impact  occurred.  This 
common  behavior  is  illustrated  in  Figure  24  together  with  a list  of 
cases  examined.  In  each  of  these  cases,  impact  with  the  floor 
occurred  at  0.66  sec  with  an  impact  velocity  of  approximately 
22  fps. 


Case  VQ  (fps)  aQ(g) 

(a)  44  3.5 

(b)  22  3.5 

(c)  22  5.0 

(d)  22  7.0 

(e)  11  5.0 


Time,  t (sec) 

FIGURE  24.  VELOCITY  PERTURBATION  - TUMBLING  PASSENGER 
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For  the  far  field  cases  where  the  impact  on  a barrier  occurs 
after  the  barrier  stops,  the  time  can  be  approximated  by  the  ratio 
of  the  barrier  distance,  x^,  to  the  initial  train  velocity,  VQ,  and 
the  impact  velocity  would  then  be  equal  to  the  initial  train 
velocity  plus  the  time  dependent  increment.  This  formulation  is 
illustrated  in  Figure  25  for  the  five  cases  examined.  The  curves 
demonstrate  the  relative  commonalty  of  the  behavior  over  a range  of 
velocity  changes  and  decelerations  of  interest. 

Additional  cases  using  the  anthropometric  categories  of 
interest  showed  that  the  motion  behavior  is  quite  insensitive  to 
the  anthropometric  categories.  Thus  a simple  impact  velocity 
formulation  can  be  generated  to  represent  all  population 
categories,  all  deceleration  levels  and  all  train  speeds  of 
interest.  A comparison  of  the  velocity  increment  variations  for 
the  anthropometric  categories  is  illustrated  in  Figure  26.  Table  4 
presents  the  floor  impact  velocities  for  these  population 
categories  and  demonstrates  the  general  insensitivity  of  this 
injury  producing  parameter  to  the  anthropometric  variances  expected 
in  intracity  trains. 

TABLE  4. -FLOOR  IMPACT  VELOCITY 


Population  Impact  Velocity 

Category*  (fps) 


Female 

Avg 

21.8 

Male 

Avg 

23.1 

Male 

Tall 

24.1 

Male 

Short 

22.0 

(23.1) 

Male 

Thin 

23.2 

— Male 

Fat 

23.0 

Child 

10  year  old 

20.1- 

Average 

Adult 

(22.4) 

Adult  except  as  noted. 
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FIGURE  25.  VELOCITY  INCREMENT  VARIATIONS  WITH  VELOCITY  AND  ACCELERATION 
FOR  THE  ADULT  FEMALE  STANDING  PASSENGER 
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FIGURE  26.  VELOCITY  INCREMENT  VARIATION  FOR  STANDING  PASSENGER 


A head  impact  velocity  of  20  fps  upon  a hard  surface,  such  as  a 
typical  train  floor,  suggests  that  some  injuries  could  be 
eliminated  or  their  severity  reduced  by  giving  special  attention  to 
the  design  of  the  floor  surface.  Figure  27  presents  .the  maximum 
acceptable  stiffness  for  an  elastically  behaving  floor- for  various 
impact  velocities  using  a severity  index  value  of  1000  sec  as  an 
arbitrary  criteria.  The  typical  impact  velocity  of  21.8  fps 
requires  a maximum  stiffness  of  approximately  1250  Ib/inch. 

4.5  Forward  Facing  Seated  Passenger 

The  term  ’’forward”  is  not  in  reference  to  the  direction  the 
train  is  traveling,  but  is  relative  to  the  direction  in  which  the 
car  velocity  decreases.  The  gross  response  of  the  forward  facing 
passenger  is  a head  forward  rotation  in  the  passenger  facing 
direction  and/or  a sliding  out  of  the  passenger's  seat. 

This  section  treats  two  limiting  models  referred  to  as  a hip 
loft  model  and  a pinned  hip  model.  It  is  expected  that  the  torso 
connection  to  the  relatively  fixed  legs  will  allow  the  hips  to  loft 
somewhat  between  the  two  extremes  allowed  by  the  two  models.  In 
this  manner  the  actual  motion  can  be  bounded.  Furthermore,  it  is 
expected  that  the  passenger  will  slide  forward  before  substantial 
torso  rotation  occurs.  This  sliding  feature  was  incorporated  into 
the  pinned  hip  model  and  the  influence  of  sliding  was  therefore 
evaluated  in  this  case.  An  examination  of  the  variability  of  the 
impact  conditions  on  an  obstacle  (such  as  a seatback  or  panel) 
immediately  forward  of  the  passenger  provided  a rationale  for  the 
specification  of  the  impact  condition  in  terms  of  the  variables  of 
interest.  This  section  is  concluded  by  a limited  examination  of 
the  influence  that  the  train  acceleration  pulse  shape  has  upon  the 
impact  velocity. 

4.5.1  Hip  Loft  Model  - The  torso  of  the  forward  facing  seated 
passenger  is  represented  by  a rigid  body  bar  whose  lower  end  (the 
hip)  is  constrained  to  move  in  the  horizontal  direction  with  the 
motion  of  the  seat  (i.e.,  the  train  car).  The  c.g.  of  the  torso  is 
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o 20  40 

Impact  Velocity,  (fps) 

FIGURE  27.  CRITICAL  STIFFNESS  FOR  FLOOR  SYSTEMS 


65 


located  a distance  D above  the  lower  end  such  that  when  the  train 
car  decelerates  the  body  bar  rotates  forward.  This  representation 
is  illustrated  in  Figure  28.  The  lower  end  of  the  body  bar  is 
defined  as  point  A,  the  c.g.  as  point  C and  the  u_pper  end  as 
point  B.  The  angle  of  rotation  9 is  positive  in  the  clockwise 
direction,  being  zero  at  the  vertical  position.  This  vertical 
position  is  the  initial  position  of  the  body  bar.  The  horizontal 
displacement  of  the  train  car  (and  hence  of  point  A)  is  described 
by  the  variable  x.  The  horizontal  displacement  of  the  c.g.  is 
described  by  the  variable  z.  The  displacement  of  the  barrier, 
which  is  originally  a distance  S from  the  point  (x  = z = 0)  is 
defined  by  the  variable  w.  The  location  and  orientation  of  the 
body  bar  at  some  preimpact  time  is  shown  by  the  bar  A'C'B*  while  its 
original  position  and  orientation  is  shown  as  ACB  in  Figure  28(b). 
The  value  of  9 = 9^  corresponds  to  the  angle  of  rotation  at  the  time 
of  impact  (see  Figure  28(a)). 

The  train  motion,  which  is  treated  as  a constant  deceleration 
motion,  is  governed  by  the  relationship 


X ( t ) = 

constant  = -a^ 

0 < t < t 

0 

— — 0 

X ( t ) = 

constant  = 0 

t < t 
€o  L 

(55) 

x(t0) 

= xo  * 

where  V , t , a and  x^  were  defined  in  Section  4.1  and  (*)  and  ( ) 
ooo  o 

represent  the  first  and  second  derivative  respectively  of  the 
variable  with  respect  to  time,  t. 

The  equations  of  motion  for  the  c.g.,  C,  are 

Wz  = -F 

and  (56) 

— 19  = M = FD  cos9 


where 

W = weight 

F = horizontal  reaction  force  at  A 
M = moment  due  to  reaction  force 
I = moment  of  inertia^ 
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FIGURE  28.  HIP  LOFT  MODEL  FOR  FORWARD  FACING  SEATED  PASSENGER 


Eliminating  the  force  F yields 

2 - - I® 

WD  cosO  \ 

The  kinematic  constraints  of  the  system  require 
z - x = D sinO 


(57) 

(58) 


which  differentiated  twice  yields 

z - x = D cosO  0 - D sinQ  (O)2^ 
The  applicable  initial  conditions  are 


2(0)  = 0 
2(0)  = Vq 
x ( 0 ) = 0 
0(0)  = 0 
0(0)  = 0. 

Combining  equations  (57)  and  (59)  results  in 


0 


cosO 


a 

+ cosO 


where 

a 


sinO  ( 0 ) 2 


x 

D 


(59) 


(60) 


It  will  be  convenient  at  this  point  to  eliminate  the  explicit 
dependence  upon  time,  t,  by  recogni2ing  that  the  angular  velocity, 
co,  is  defined  as 


co 


dO 

dt 


so  that 
0 = 


dco  dco 

= dt  = d0 


dO  (1  dco 

dt  dO 


1/2 


d(co2)  . 
dO 


(61) 


(62) 


Further  simplifications  will  result  if  the  angular  velocity  is  made 
dimensionless  as 


(63) 
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and  introducing  the  following  parameter 


K 


(64) 


Equation  (60)  can  be  rewritten  as 


X cosQ  + - sinQ  X2  = k , 


0 < t < t 


o 


(65) 


and 


(66) 


for  each  of  the  two  constant  acceleration  phases  of  the  train 
motion . 

A limiting  solution  corresponding  to  an  instantaneous  stopping 
of  the  train  car  ( < -*■  «>)  results  in  the  replacement  of  equation 
(65)  by  a new  boundary  condition  on  X(o),  namely  X(o)  = 1,  that  is 
the  translational  momentum  is  converted  instantaneously  to  rotation 
momentum  as  Doo  = V . Thus  equation  (66)  subject  to  the  boundary 
condition  can  be  integrated  to  yield 


A nominal  value  of  9^  (the  angle  at  which  impact  occurs)  of 
45  degrees  was  selected  as  being  representative  of  the  forward 
facing  seated  passenger.  Two  values  of  a were  selected  as  repre- 
senting the  range  of  interest,  these  were  a = 0.126  and  a = 0.3* 
The  latter  value  was  considered  to  be  the  most  appropriate.  The 
angular  velocity  variation  with  the  angle  of  rotation  are  shown  in 
Figures  29  and  30  and  indicate  that  for  the  limiting  case  as  k + «> 
the  influence  of  the  parameter  a is  small. 

The  time  of  impact,  or  the  time  details  of  this  passenger 
motion  will  be  of  some  value,  hence  it  will  be  convenient  to 
introduce  a dimensionless  time  variable  as 


cos  9 + a 


(67) 


tv 


o 


(68) 


T 


D 
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FIGURE  29-  ANGULAR  VELOCITY  - ROTATION  RELATIONSHIP 
FOR  HIP  LOFT  MODEL,  a = 0.126 
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FIGURE  30.  ANGULAR  VELOCITY  - ROTATION  RELATIONSHIP 
FOR  HIP  LOFT  MODEL  - a = 0.3  - 
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Thus  the  dimensionless  form  of  equation  (61)  together  with 
equation  (67)  yields 
Q 


■l 


cos  Q + a 
1 + a 


dO. 


(69) 


For  the  general  case  ( x < «)  the  deceleration  phase  of  the 
train  ends  when 


t 

o o 

To  = D 

and  using  the  results  of  Section  4.1  yields 
V2 

T, 


(70) 


_o ]_ 

lo  aQD  = < 


(71) 


In  general  the  time  is  given  as 
d 

T 


■I 


dQ 

TToT 


(72) 


Equation  (65)  was  solved  numerically  (subject  to  the  boundary 
condition  A(0)=0)  together  with  equation  (72)  until  t = t . At 
that  time  point  A = A ^ and  9 = Q ^ and  the  solution  for  the  second 
phase  becomes,  from  equation  (66) 


and 


cos  9,  + a 


A = A. 


cos^S*  a 


(73) 


T = To  + 


J 

/ 


cos  9 + a 
2 

cos  0 i + a 


dQ 


(74) 
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The  time  history  was  obtained  by  numerical  integration  of  equa- 
tion (74).  The  results  of  these  integrations  are  presented  in 
Figures  29  and  30  for  the  angular  velocity  and  in  Figure  31  for  the 
time  details.  The  time  details  are  only  presented  for  a = 0.126. 
Four  representative  values  of  < of  0.25,  0.50,  0.75  and  1.00  were 
s-elected . 

The  above  solutions  define  the  motion  history  of  the  body  bar 
and  can  now  be  used  to  evaluate  the  velocity  of  the  potential  impact 
point  (point  B of  Figure  28).  The  dimensionless  velocity  of  the 
train  case  is 


o 


and 

X = 0 


> 


0 < X < T 

- - 0 


T < T 

0 — o 


Whereas  the  dimensionless  velocity  of  the  c.g 
u 4 

= v7 


is 


(75) 


(76) 


and  the  relative  velocity  between  the  c.g.  and  the  barrier  is  from 
the  derivative  of  equation  (58) 

e>zr  = “ x = XcosQ.  (7? ) 

The  relative  velocity  dependence  on  the  angle  of  rotation  is 
presented  in  Figure  32.  The  time  details  of  the  various  velocities 
are  presented  (for  a = 0.126)  in  Figure  33*  It  is  clear  from  these 
results  that  the  horizontal  velocity  of  the  c.g.  does  not  change 
significantly,  such  that  the  relative  velocity  is  controlled 
primarily  by  the  motion  of  the  train  car.  If  the  body  component 
impact  occurs  during  the  second  phase  (x  < t ) of  the  motion,  and 
it  frequently  can  and  does  (see  Figure  32  for  9 = 9^),  the  relative 
velocities  are  high. 

The  details  of  an  impact  upon  an  obstacle  located  immediately 
forward  of  the  passenger  is  illustrated  in  Figure  347a).  Since 
there  is  no  vertical  force  acting  upon  the  c.g.  it  moves  only  in  the 
horizontal  direction. 
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FIGURE  31.  ROTATION  HISTORY  FOR  HIP  LOFT  MODEL 
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FIGURE  32.  CENTER  OF  GRAVITY  VELOCITY  FOR  HIP  LOFT  MODEL 
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FIGURE  33.  CENTER  OF  GRAVITY  VELOCITY  HISTORY  FOR  HIP  LOFT  MODEL 
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FIGURE  34.  IMPACT  DETAILS  FOR  HIP  LOFT  MODEL 


Thus  the  hip  (point  A)  moves  upward  (to  position  A')  during  the 
forward  rotating  motion.  The  velocity  components  of  the  body  bar 
at  the  angle  of  impact  (0^)  are  illustrated  in  Figure  3Mb)  rela- 
tive to  the  car  interior.  The  velocity  of  the  c.g.  is  V . The 
velocity  of  point  A'  is  composed  of  two  components;  the  c.g. 
component  and  the  rotational  component  where  = wD. 


The  resultant  velocity  of  these  two  components  is  VA  which  is 
vertical  in  accordance  to  the  applied  constraint  on  the  lower  end 
of  the  body  bar.  The  velocity  component  of  point  Bf  (the  impact 
point)  is  also  composed  of  a c.g.  contribution  and  a rotational 
contribution,  Vg  where  Vg  = ojD.  From  the  definition  of  the 

dimension  H in  Figure  28(b)  it  is  clear  that  H = D + D’. 

* 

The  resultant  velocity  of  point  BT  is  Vg.  The  normal  impact 
velocity  will  depend  upon  the  orientation  of  the  impact  surface. 

This  feature  is  illustrated  in  Figure  34(c)  for  three  surfaces 

# 

normal  respectively  to  the  velocity  vectors  Vzr,  Vg  and  Vg.  The 
corresponding  resultant  impact  velocities,  V^ , V^  and  Vg  are  all 
about  equal  in  magnitude  hence  the  well  defined  value  Vg  can  be 

used  and  the  orientation  of  the  surface  need  not  be  considered 

# 

explicitly.  The  value  of  Vg  in  dimensionless  terms  is 


i . 

imax 


zr 


D T 

+ Ag  cosO^ 


D ’ 

Ag  sinGg 


(78) 


The  dependence  of  the  contributing  parameters  <t>  and  A are 

i 

presented  in  Figure  35  for  the  case  D'  = D as  a function  of  k.  The 
two  distinct  regions  correspond  to  an  impact  occurring  in  each  of 
the  two  phases  of  the  railcars  motion.  Furthermore,  the  parameter 
k can  be  redefined  using  the  results  of  Section  4.1  as 


(79) 


o * 


and  thus  is  related  simply  to  the  car  stopping  distance  x 


D=D 


FIGURE  35.  INFLUENCE  OF  ACCELERATION  ON  VELOCITY  COMPONENTS 
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4.5.2  Pinned  Hip  Model  - The  pinned  hip  model  restrains  the 
lower  end  of  the  body  bar  (point  A)  from  moving  upward  during  the 
rotation  of  the  body  bar.  Thus  an  additional  equation  of  motion  for 
the  vertical  movement  of  c.g.  is  required,  namely 

Wy  = G (80) 

where 

y = the  downward  displacement  of  the  c.g. 

G = the  vertical  reaction  force  at  A. 

The  moment  equation  (in  equation  56)  must  be  generalized  and 
rewritten  as 


M = G(z-x)  - yF 

and  the  kinematic  constraints  must  be  expanded  to  include 
y = D cosQ. 


(81) 

(82) 


Following  a similar  procedure  as  before  yields  the  basic  differen- 
tial equation 


0 ( 1 + a)  = ” § cosQ*  ( 83 ) 

When  this  equation  is  expressed  in  terms  of  the  dimensionless 
variables  for  each  of  the  two  phases  of  the  train  car  motion,  these 
results  are  obtained 


A ^ (1  + a)  = k cosO,  0 < T 1 T0  (84) 

and 

X ( 1 + a)  = 0 , T0  <_  x (85) 

or 

x = 0 > To  i T-  (86) 

If  the  limiting  solution  for  k is  sought  the  initial  condi- 
tion is  as  before  A(0)  = 1 and  the  differential  equation  (equation 
(86))  for  tq  <_  t applies.  Integration  yields 

A constant  = 1 , 0 £ t 

and  from  equation  (72)  further  integration  yields 
0 = x. 
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The  more  general  case  (<  < <»)  can  also  be  integrated  directly 
recalling  that  the  boundary  condition  A(0)  = 0 applies  for  equa- 
tion ( 84 ) . Thus 

A = A*  V sinO  , 0 _<  t _<  tq  ‘ (87) 

where 

A*  = 


1 + a 


(88) 


The  corresponding  results  for  this  phase  are  presented  in 
Figure  36.  The  time  details  are  obtained  numerically  from 


'■irf  — 


V sinO 


0 < t < x 


(89) 


however,  auxiliary  solution  can  be  utilized,  namely,  let 


FIGURE  36.  ANGULAR  VELOCITY  - ROTATION  RELATIONSHIP 
FOR  PINNED  HIP  MODEL 
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This  function  was  evaluated  and  the  results  are  presented  in  Fig- 
ure 37.  In  particular  when  t = tq  and  Q(t  ) = 0^ 

T(  0 1 ) = A*T0  = -f-  (91) 

the  value  of  = X(O-j)  is  then  given  as 

X-|  = X*  VsinO  -|  • (92) 


The  solution  for  the  second  phase  xQ  < t is  simply 
X = constant  = X^ 


and 


T = To  + 


( ©-0 1 ) 


To  -6  1 


(93) 


The  complete  angular  velocity  and  time  details  for  this  case  are 
presented  in  Figures  36  and  38. 


The  dimensionless  relative  horizontal  velocity,  6 is 

zr 

identical  (see  equation  (77))  to  that  for  the  hip  loft  model. 
However,  for  this  model  the  vertical  velocity,  must  also  be 

defined,  namely 

6 - t-  = • X sin©  (94) 

y o 

These  velocity  components  are  presented  in  Figure  39. 

4.5*3  Influence  of  Slip  - The  likelihood  that  a forward  facing 
seated  passenger  will  slide  forward  in  his  or  her  seat  before 
substantial  torso  rotation  will  occur,  appears  to  be  great.  Thus 
this  feature  was  incorporated  into  the  pinned  hip  model.  It  was 
hoped  that  such  an  evaluation  would  yield  a meaningful  measure  of 
this  effect.  The  model  assumes  that  the  passenger's  horizontal 
motion  remains  constant  at  its  initial  value,  V , after  the  train 
car  deceleration  starts  (t  = 0)  and  until  a relative  distance,  s, 
is  traversed  at  t = t'.  This  distance  represents  the  knee-to-seat- 
back  clearance  which  exists  and  thus  it  becomes  a parameter  in  the 
problem.  During  this  period  the  passenger  moves  a distance  zf  and 
the  train  moves  a distance  x'. 
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FIGURE  38.  ROTATION  HISTORY  FOR  PINNED  HIP  MODEL 
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FIGURE  39.  CENTER  OF  GRAVITY  VELOCITY  COMPONENTS  FOR  PINNED  HIP  MODEL 


These  distances  are  obtained  from  the  respective  acceleration 
histories  as 

(95) 


Z'  = V0t 

and  ant2 

x'  = VQt  - 4-. 


But  at  t = t ' 

s = z T - x ’ 

hence 


2s 

a 


(96) 

(97) 

(98) 


or  in  dimensionless  terms 

■ tf  ■ V¥" 


(99) 


where 

a = §’  (100) 

When  knee  impact  occurs  the  model  assumes  that  the  linear 
momentum  is  converted  instantaneously  to  angular  momentum  and 
rotation  of  the  vertical  body  bar  begins.  Thus  the  new  initial 
conditions  for  t = t*  are 


0 ( t ’ ) = 0 


0) 


(t»)  = 


AV 

D 


( 101 ) 


where  aV  is  the  instantaneous  velocity  change  that  point  A experi- 
ences at  the  time  of  impact  (see  Figure  40(a)).  This  velocity 
change  is  simply  aQt ’ provided  that  tr  _<  tQ.  This  time  limitation 
is  usually  satisfied  since  s will  generally  be  much  less  than  x . 
In  dimensionless  terms  the  second  boundary  condition  becomes 


= T*  K a 


(102) 


Before  the  solution  is  presented  it  will  be  appropriate  to 
recognize  that  due  to  the  forward  slippage  of  the,  point  A the 
clearance  between  the  passenger  and  the  barrier  which  was 
originally  S has  been  reduced  to  a value  S - s.  Thus  the  amount  of 
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FIGURE  40.  HIP  SLIP  CHARACTERISTICS  FOR  PINNED  HIP  MODEL 
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rotation  which  would  have  occurred  without  slippage,  9!,  has  been 
reduced  to  9!  with  slippage.  This  effect  is  illustrated  in  Fig- 
ure 40(b).  This  angle  is  given  as 

S-s 


9?  = arc  sin 


D+D 


whereas  without  slippage 

S 


9^  = arc  sin 


D+D 


(103) 


(104) 


Using  the  conditions  D*  = D and  9^  = 45  degrees  requires  that 
S = /2D,  thus 


9j  = arc  sin 


j ( /~2-'cT) 


(105) 


Since  the  value  of  D is  approximately  1 ft  values  of  a = 0.25  and 
0.5  were  selected  as  being  appropriate  for  this  evaluation.  These 
correspond  to  values  of  9J  Of  38.4  degrees  and  27.2  degrees  respec- 
tively (or  0.621  and  0.475  rad). 


The  basic  equation  (equation  84)  can  now  be  integrated  to 


yield 


2 2 
* = 4 


2 k sinQ 
( 1 + a) 


which  with  the  use  of  equations  (99)  and  (102)  becomes 

A = ~\J  2 k ~\J"° 


’G  + 


sinQ 
1 +a 


T 1 _<  T _<  XQ » 


Define,  as  before 
xi  = 


then  for  the  second  phase 
A = constant  = A-| 


To  1 T* 


(106) 


(107) 


(108) 


(109) 


A representative  solution  is  illustrated  in  Figure  41.  The  time 
details  must  be  evaluated  numerically  for  both  phases  from  the 
basic  definition,  viz« 


T = T 


9 

■*/■ 


d0 

aT9T 


( 110) 
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9 (rad) 

FIGURE  41.  INFLUENCE  OF  SLIP  ON  ANGULAR  VELOCITY  - 
PINNED  HIP  MODEL 


The  velocity  components  are  given  by  the  same  basic  definitions  as 
before,  however  their  values  are  different.  These  results  are 
given  in  Figure  42  for  the  representative  solution. 

4.5.4  Variability  of  Impact  Conditions  - The  previous  three 
subsections  have  developed  results  with  which  to  evaluate  the 
impact  conditions  for  the  hip  loft  model,  the  pinned  hip  model  and 
t-o  measure  the  influence  of  slip  for  the  pinned  hip  model.  . For  the 
nonslip  condition  the  body  component  impact  environment  should  be 
bounded  by  the  two  models.  The  amount  of  slip  which  will  occur  will 
be  less  than  that  permitted  by  the  above  slip  conditions  used  in  the 
slip  version  of  the  pinned  hip  model.  The  value  of  a will  vary  with 
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FIGURE  H2.  INFLUENCE  OF  SLIP  ON  CENTER  OF  GRAVITY 
VELOCITY  COMPONENTS  - PINNED  HIP  MODEL 


both  car  seat  design  and  spacing  and  anthropometric  variances.  It 
is  hoped  that  the  two  models  in  addition  to  the  range  of  a selected, 
will  represent  a narrow  enough  band  of  impact  conditions  such  that 
this  collection  of  information  can  be  used  to  est-ablish  some 
nominal  impact  environment. 

Figure  43(a)  presents,  for  a given  set  of  conditions,  the 

velocity  vectors  corresponding  to  the  hip  loft  model,  where  d. 

imax 

was  defined  in  Section  4.5.1.  Figure  43(b)  presents  the  velocity 
vectors  corresponding  to  the  pinned  hip  model  without  slip  ( a = 0) 
and  the  resultant  vectors  for  different  slip  distances,  a.  The 
vector  a is  identical  to  the  resultant  vector  for  the  corresponding 
hip  loft  model.  The  magnitude  of  the  impact  vector  for  the  two 
nonslip  cases  are  approximately  equal  and  not  much  smaller  than  the 
magnitudes  of  those  vectors  for  the  various  slip  distances.  It 
would  appear  that  both  the  magnitudes  and  the  direction  are  similar 
enough  to  make  use  of  some  averaged  condition  and  to  threat  the 
calculated  variance  as  just  another  variance  contribution  is  a 
situation  where  many  variances  will  exist.  Figures  44  amd  45 
present  the  maximum  impact  velocity  dependence  on  the  parameter  < 
for  the  cases  examined  above  and  for  one  set  of  anthropometric  and 
seat  spacing  values. 

4.5.5  Influence  of  Pulse  Shape  - The  pinned  hip  model  (without 
slip)  was  used  to  evaluate  the  influence  that  train  car  decelera- 
tion pulse  shape  has  upon  the  body  component  impact  environment  for 
the  forward  facing  seated  passenger.  This  model  was  -selected 
because  of  its  relative  simplicity,  however  no  specific  derivations 
are  presented.  The  pulse  shape  that  was  selected  was  the  double 
step  pulse  (see  Section  4.1)  with  the  parameter  values  6 = 0.5  and 
t<|  = 0.5;  thus  g = 0.417  and  y = 4/3  (see  Figure  14).  Equivalency 
in  this  evaluation  was  based  upon  conserving  both  the  change  in 
velocity,  V , and  the  stopping  distance,  x . The  conservation  of 
the  stopping  distance  was  selected  in  view  of  the  importance  that 
the  parameter  k plays  in  defining  the  severity  of  this  impact 
environment . 
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s = VFd 


FIGURE  43.  VECTOR  DIAGRAMS  FOR  IMPACT  VELOCITIES 
OF  SEATED  PASSENGER  FACING  FORWARD 
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FIGURE  44.  IMPACT  VELOCITY  DEPENDENCE  ON  MODEL 
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FIGURE  45.  IMPACT  VELOCITY  DEPENDENCE  ON  SLIP 
PINNED  HIP  MODEL 


94 


This  equivalency  requires  that  Bt  ( = x /V  ) remain  constant  for 

ooo 

the  two  pulses,  thus  the  pulse  duration  for  the  double  step  pulse 
must  be  1.2  times  larger  than  the  reference  step  pulse  duration  and 
the  average  acceleration  must  be  1.2  times  less  than  that  of  the 
reference  step  pulse. 

The  results  for  the  angular  velocity  for  these  two  equivalent 

cases  are  presented  in  Figure  46  and  indicate  that  in  general  the 

differences  are  not  large,  in  particular  if  6^  = 0.8  rad.  The 

influence  of  the  pulse  shape  is  more  properly  evaluated  by 

examining  their  influence  on  the  impact  velocity  d • . This 

imax 

comparison  is  presented  in  Figure  47  and  clearly  demonstrates  that 
the  influence  of  the  pulse  shape  is  quite  small  and  generally  less 
than  the  variations  or  uncertainties  introduced  by  the  degree  of 
hip  loft  or  slip,  or  by  anthropometric  uncertainties. 

4.6  Backward  Facing  Seated  Passenger 

The  backward  facing  seated  passenger  will  be  pressed  back  into 
his  or  her  seat  due  to  the  sudden  change  in  motion  brought  about  by 
the  railcar  crash.  Therefore,  the  motion  of  passengers  in  this 
passenger  configuration  will  be  generally  restrained  and  no  signif- 
icant impact  events  are  expected.  The  primary,  potentially 
detrimental,  response  will  be  the  so-called  whiplash  or  neck 
flexion  response. 

Flexion  response  data^  for  the  50th  percentile  adult  male  is 
presented  in  Figure  48  and  represents  the  basis  for  the  following 
analysis.  A pinned  bar  model  was  used  for  this  analysis  where  the 
bar  represents  the  head  and  the  lower  pinned  end  (point  A in  the 
previous  analysis)  represents  the  neck  connection.  The  neck  will 
offer  a resistance  to  rotation  as  illustrated  in  Figure  48  and  the 
onset  of  injury  will  occur  after  the  head  has  rotated  approximately 
70  to  80  degrees.  A linear  rotational  resistance,  B,  was  assumed 
for  this  model  and  a nominal  value  of  50  ft  lb/rad  was  selected. 


"Biomechanics  and  Its  Application  to  Automotive  Design",  Society 
of  Automotive  Engineers,  p 49,  1973* 
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FIGURE  46.  INFLUENCE  OF  PULSE  SHAPE  UPON  ANGUL-AR  VELOCITY 
PINNED  HIP  MODEL  FOR  FORWARD  FACING  SEATED 
PASSENGER 
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FIGURE  47.  INFLUENCE  OF  PULSE  SHAPE  ON  IMPACT  VELOCITY 
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Flexion  Moment  about  Occipital  Condyles  (ft-lb) 


FIGURE  48.  FLEXION  RESPONSE  ENVELOPES 

FOR  50TH  PERCENTILE  ADULT  MALE 
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No  attempt  was  made  to  simulate  the  increased  resistance  which 
occurs  at  large  rotations  since  the  energy  involved  in  this  addi- 
tional rotation  is  not  excessive  and  is  small  compared  to  the 
uncertainty  of  the  rotation  resistance  at  smaller  angles  of 
rotation.  Furthermore,  a conservative  critical  angle,  Qc,  of 
65  degrees  (1.135  rad)  was  selected  as  the  injury  threshold.  In 
this  manner  the  analysis  is  rather  conservative.  The  other 
pertinent  anthropometric  data  were  presented  in  Section  4.3. 

The  basic  equation  and  initial  conditions  used  for  the  pinned 
hip  model  (see  Section  4.5.2)  apply  for  this  model,  except  that  the 
moment  equation  must  be  expanded  to  include  the  rotational  resis- 
tance term,  namely 

M = G(z-x)  - yF  - BO  (111) 


and  the  basic  differential  equation  becomes 


Q(1+a)  = - ^ cosQ  - 

U WD 


( 112) 


which  when  rewritten  in  its  dimensionless  form  for  the  two  train 
motion  phases  become 


X (1+o0  = kcosQ  — 2ke0  , 


and 


X ( 1 + a)  = 2 keO 


0 _<  T _<  T0 


To  - T 


where 


e = 


B 


(113) 

(114) 

(115) 


2WV*~k  2a^  WD 

o o 


Integration  of  equation  (113)  yields 

x = x/f^Vsine  • e°2  • 


( 116) 


Note  that  the  parameter  e involves  the  acceleration  aQ.  If  k is 
small,  which  will  be  the  case  for  long  duration  acceleration  pulses 
(t  large),  the  maximum  neck  rotation  will  occur  during  the  first 
phase.  The  maximum  value  of  0 can  then  be  obtained  simply  from 


e 


sinO 


(117) 
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or  more  importantly  a critical  value  of  e,  say  £q  can  be  evaluated 
corresponding  to  the  critical  (threshold)  value  of  ©,  namely  © . 
In  this  manner  a critical  or  limiting  value  of  aQ  can  be  determined, 
namely 


B 


2ecWD 


(118) 


For  the  anthropometric  values  used  the  limiting  acceleration  is 
11.8  g and  ec  = 0.705.  Since  most  train  crash  environments  are 
expected  to  have  average  decelerations  considerably  less  than  this 
amount  it  would  appear  the  whiplash  injuries  will  be  minimal. 


For  extremely  high  train  decelerations  (i.e.,  the  instantane- 
ous stopping)  the  limiting  solution  k -*  °°  will  be  of  value.  As 
before  the  second  phase  equation  (t  £ x)  applies  with  the  special 
boundary  condition  (X(0)  = 1) 


A =~\/l  - 62  (119) 

V 1 + a 

where 

y = 2 k e = — ^ ' (120) 

WV^ 

o 


There  exists  a corresponding  critical  value  for  velocity  change 
VQ  Q which  will  limit  © to  the  critical  value  ©c,  namely 


V 


ecB 


( 1+  cOW 


( 121  ) 


Since  a = 0.3*  Yc  = 1.01  and  VQ  = 12.6  fps. 


Using  the  general  relationship  (equation  (72))  for  the  dimen- 
sionless time,  we  can  again  define 


A ( T0)  = A1 


( 122) 


and 


© ( A1 ) = ©1 


( 123) 
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and  integrate  the  basic  equation  for  the  second  phase, 
results  in 


This 


t < i, 

o - 


(124) 


Typical  results  of  the  angular  history  of  the  head  are  presented  in 
Figure  49.  Maximum  rotation  contours  are  presented  in  Figure  50(a) 
and  when  compared  with  Qc  permit  the  mapping  of  "safe"  and  "unsafe" 
domains  in  the  k - e plane.  These  domains  are  presented  in 
Figure  50(b).  The  safe  region,  when  mapped  in  the  train  crash 
domain  (see  Figure  51)  provides  an  immediate  impression  of  the  low 
hazard  that  backward  facing  seated  passengers  are  exposed  to  in 
typical  crash  environments. 


4.7  Sideward  Facing  Seated  Passenger 

The  last  of  the  passenger  configurations  studied  is  that  of 
the  sideward  facing  seated  passenger.  The  body  bar  representation 
was  again  selected  to  describe  the  primary  motion  characteristics 
of  this  passenger  configuration  to  the  sudden  deceleration  of  the 
train  car.  Figure  52(a)  illustrates  the  passenger  motion.  The 
whole  body  c.g.  for  this  seated  condition  is  just  above  the 
apparent  point  of  rotation  and  since  no  lofting  effect  is  expected 
a pinned  body  bar  (at  point  A)  representation  was  used.  The  basic 
equations  are  identical  to  those  presented  in  Section  4.5.2  for  the 
forward  seated  passenger.  However,  in  this  configuration  the 
anthropometric  variables  are  numerically  different  and  thus  the 
corresponding  parameter  values  are  substantially  different.  The 
moment  of  inertia  of  the  seated  whole  body  is  larger  and  the  charac- 
teristic dimension  D is  somewhat  smaller;  thus  a = 1 1 to  17  and  the 
range  of  k of  interest  is  approximately  0.02  to  0.10.  Furthermore 
the  H/D  ratio  is  larger,  approximately  10. 

The  dimensionless  angular  velocity  is  smaller  due  to  the 
smaller  value  of  D.  Typical  results  for  the  angular  velocity  are 
presented  in  Figure  53.  The  dimensionless  impact  velocity  ^max  is 
approximately  the  same  as  for  the  forward  facing  seated  passenger. 
The  results  are  presented  in  Figure  54  and  a range  of  a values  and 
impact  rotation  angles  of  0^  = 0.4,  0.5,  and  0.6  rad. 
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FIGURE  49.  ANGULAR  VELOCITY  HISTORY  FOR  NECK  MODEL 
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FIGURE  50.  RESPONSE  CHARACTERISTICS  OF  NECK  MODEL 
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FIGURE  51.  SAFE  REGION  IN  TRAIN  CRASH  DOMAIN 

FOR  BACKWARD  FACING  SEATED  PASSENGERS 
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FIGURE  52.  PINNED  HIP  MODEL  FOR  SIDEWARD  FACING  SEATED  PASSENGER 
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FIGURE  53.  ANGULAR  VELOCITY  FOR  SIDEWARD  SEATED  PASSENGER 
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FIGURE 


54.  IMPACT  VELOCITY  FOR  SIDEWARD  SEATED  PASSENGER 
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4.8  Multiple  Passenger  Interactions 

It  is  quite  evident  that  whenever  train  cars  are  crowded,  as 
they  frequently  are  during  the  peak  hours,  that  standing  passengers 
in  particular  will  bump  into  each  other  during  a crash  exposure. 
Whereas  an  isolated  passenger  may  impact  a barrier  once, 
interactions  with  nearby  passengers  may  cause  additional  impacts, 
involving  the  same  body  component.  The  secondary  passengers  who 
are  isolated  from  the  barrier  (wall  or  floor)  by  another  passenger 
may  only  experience  a crush  environment  and  may  not  be  subjected  to 
any  significant  impact  condition.  The  concern  then  is  for  the 
primary  passenger  who  is  adjacent  to  the  barrier  and  who  may 
experience  multiple  impacts. 

The  approach  used  herein  to  attempt  to  obtain  a measure  of  the 
increased  hazard  that  passenger  interactions  may  cause  to  a given 
passenger  is  to  examine  the  multiple  impact  conditions  which  may 
result  from  a short  series  of  multiple  passenger  interactions  and 
to  examine  the  potential  increase  that  such  impact  environments  may 
have  upon  injury  criteria,  such  as  the  severity  index.  It  is  clear, 
at  the  onset,  that  the  multiple  impact  scenario  selected  is  but  one 
of  many  and  that  it  does  not  necessarily  represent  a very  severe 
scenario.  It  is  recognized  that  these  multiple  impacts  will  occur 
over  a time  interval  for  which  no  injury  criteria  are  directly 
applicable.  When  a so-called  total  severity  index  is  formulated  by 
summing  the  contribution  from  individual  impacts,  it  is  understood 
that  any  interpretation  of  the  final  result  can  only  be  used  as  a 
very  rough  indication  of  the  effect  of  multiple  impacts.  Two 
barrier  types  (the  perfect  absorber  and  the  elastic  response 
barrier)  were  selected  for  analysis  and  the  passenger-passenger 
interaction  was  assumed  to  be  inelastic  such  that  the  masses  asso- 
ciated with  the  passengers  remained  in  contact.  The  total  momentum 
was  conserved. 

The  path  diagram  for  the  first  case,  the  perfect  absorber) is 
illustrated  in  Figure  55(A).  The  primary  passenger,  PI  of  weight 
W i = W impacts  the  barrier  at  a velocity  VQ.  Thus  the  severity 
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index  takes  on  a value  (SI)  which  is  proportional  to  V.  = V and 

° 1 c 10 

inversely  proportional  to  W . The  remaining  parameters  of  the 
severity  index  are  assumed  to  be  constant.  The  body  component  is 
assumed  to  separate  slightly  from  the  barrier  and  come  to  rest.  The 
second  passenger  P2  with  the  weight  and  velocity  indicated,  impacts 
inelas ti cally  upon  the  primary  passenger  and  both  move  at  the 
velocity  1/2  V , the  impact  velocity  for  the  second  impact.  Using 
the  increased  weight  and  the  decreased  velocity,  yields  a second 
severity  indxex  of  0 . 1 8 (SI)  such  that  current  total  severity 
index  is  1.18  (SI)  . If  only  the  original  body  component  weight  W 
is  involved  (such  as  it  might  be  for  the  head),  the  current  total 
severity  index  is  increased  at  this  time  to  1.5  (SI)  . The  former 
value  appears  small  enough  to  indicate  that  little  additional 
hazard  is  produced  by  multiple  passenger  interactions  (note  that 
five  interactions  yield  a total  severity  index  of  only  1.29  (SI)  . 
The  latter  total  severity  index  (i.e.,  the  head  impact)  value  is 
sufficiently  large  to  be  of  some  concern.  While  it  may  be  unlikely 
that  more  than  two  head  impacts  may  occur  at  one  general  location  of 
the  head,  it  is  difficult  to  rule  out  a double  head  impact  event. 

A similar  evaluation  was  made  using  an  elastic  responding 
surface  (see  Figure  55(b)).  The  body  component  leaves  the  barrier 
with  a velocity  VQ  after  the  first  impact  such  that  the  first 
passenger-passenger  impact  results  in  both  passengers  being  at  rest 
a short  distance  from  the  barrier.  The  collision  of  a third  passen- 
ger then  causes  a barrier  impact  at  a velocity  of  1/3  V . The 
current  total  severity  index  only  increases  to  1.06  (SI)  for  the 
same  body  component  weight.  If  a larger  weight  were  used  such  as  3W 
then  the  total  severity  index  would  be  even  smaller.  In  this  second 
case  the  passenger  interaction  effect  is  clearly  not  significant, 
at  least  not  with  respect  to  injury  due  to  body  component  impact 
phenomena . 


4.9  Nomenclature 

a acceleration 

aQ  average  acceleration 

* 

aQ  peak  acceleration 
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a1 

specific  value  of  acceleration  for  double  step  pulse 

a ’ 

retarding  acceleration 

ao,c 

critical  value  of  acceleration  for  flexion  response 

B 

rotational  resistance  for  head  rotation 

D,Df 

characteristic  body  dimensions 

F 

horizontal  force 

G 

vertical  force 

H 

body  dimension,  height  of  standing  passenger 

H 

mean  height  of  standing  passenger 

I 

moment  of  inertia 

I 

mean  value  of  moment  of  inertia  for  standing  passenger 

k 

spring  constant 

L 

body  dimension 

r 

mean  value  of  body  dimension  for  standing  passenger 

M 

moment 

/ 

s 

slip  distance 

S 

body  dimension,  obstacle  spacing 

S 

mean  value  of  body  dimension  for  standing  passenger 

SI , ( DI ) 0 

severity  index 

t 

time 

pulse  duration 

t1,t2 

time  parameters  for  pulse  shape 

t ’ 

time  at  which  slip  ceases 

t . 
1 

time  of  impact 

V 

velocity 

V 

0 

change  in  velocity 

V. 

l 

impact  velocity 

< 

o 

o 

critical  value  of  velocity  for  flexion  response 

Ill 


AV 

change  in  velocity  when  slip  ceases 

4Vx 

velocity  perturbation 

VA’VB 

* * 

VVB 

velocity  contribution  due  to  rotation 
total  velocity 

VVt 

specific  velocity  components 

Vzr 

relative  velocity  in  horizontal  direction 

w 

body  dimension,  displacement 

w 

mean  value  of  body  dimension  for  standing  passenger 

W 

whole  body  or  body  component  weight 

W 

mean  value  of  whole  body  weight  for  standing  passenger 

w1,w2,w3 

weight 

WT 

total  weight 

X 

displacement 

xo 

stopping  distance 

xb 

* 

X 

barrier  distance 
passenger  location 

*3’*c.g. 

velocity  components  for  tumbling  block 

x» 

train  displacement  when  slip  ceases 

y 

displacement 

z 

displacement 

z ’ 

passenger  displacement  when  slip  ceases 

a 

inertia  parameter,  acceleration  parameter 

a 

c 

specific  acceleration  parameter 

6 

stopping  distance  parameter 

6 

acceleration  parameter 

y 

acceleration  parameter,  parameter  for  flexion  response 
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^c 

£ 

£ 

C 

a 

T 

T 

O 
T ' 

T1  ’ T2 
X 

X' 

X 

o 

X* 


0) 

Q 


Q . , Q ! 
1 ’ l 


C 

* 

5 


X 


"i 

4. 

imax 


0 . d ,<& 
z ’ zr’  y 


critical  value  of  parameter 

acceleration  parameter  for  flexion  response 

critical  value 

scaling  factor,  spacing  parameter 
time 

pulse  duration 

time  when  slip  ceases 

time  parameter  for  pulse  shape 

scaling  parameter,  angular  velocity 

specific  value 

specific  value  when  slip  ceases 

parameter 

angular  velocity 

rotation 

rotation  at  impact 
specific  value 
critical  value 
displacement 
passenger  location 
velocity  of  train 
impact  velocity 
maximum  impact  velocity 
velocity  components 
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5.  PASSENGER  INJURY  PREDICTION  METHODOLOGY 


5.1  General  Approach 

A computer  program  based  upon  the  preceding  motion  and 
response  analyses  was  prepared  to  conveniently  evaluate  injury  pro- 
files of  passenger  groups  experiencing  an  intracity  train  crash. 
The  objective  of  the  methodology  is  to  evaluate  the  potential  con- 
sequences of  such  a crash  as  a function  of  the  crash  environment, 
car  design,  car  interior  component  design  and  arrangements,  and 
other  factors.  The  computer  program  "Intra-City  Crash  Injury 
Analysis"  is  presented  and  documented  in  the  appendix  of  this 
volume.  The  documentation  includes  a description  of  the  required 
input  data  format  as  well  as  the  output  data  format.  This  section 
presents  the  overall  structure  of  the  code  and  the  general  philos- 
ophy used  in  treating  each  of  the  four  basic  passenger  configura- 
tions. The  results  of  a sample  problem  are  then  presented  at  the 
end  of  the  section. 

The  basic  structure  of  the  passenger  injury  prediction 
methodology  is  presented  in  Figure  56  in  terms  of  11  activities. 
The  activities  define,  in  general,  the  specific  tasks  to  be 
performed,  data  which  are  to  be  provided,  or  conclusions  reached. 
This  structure  also  illustrates  the  general  flow  of  the  calculation 
procedure  and  the  interrelationships  between  certain  activities. 

Activity  1 is  a statement  of  the  passenger  distribution  with 
respect  to  the  number  of  passengers.  This  information  is  given  for 
each  kth  passenger  configuration  (K  = 1,4),  in  each  jth  car  in  each 
ith  consist.  One  or  two  consists  can  be  specified  defining  a 
consist-consist  crash  or  an  abutment  (or  other  nonpassenger 
entity )-consist  crash  situation. 

Activities  2 and  3 deal  with  an  analysis  of  the  gross  motion 
(i.e.,  car  frame  motion)  of  the  individual  cars  of  the  consists  and 
provides  the  basic  excitation  used  for  the  passenger  motion 
analyses.  These  data  include  the  absolute  values  of  both  the 
change  in  velocity  and  the  average  acceleration  for  each 
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FIGURE  56.  PASSENGER  INJURY  PREDICTION  METHODOLOGY 


car  for  each  consist.  These  values  correspond  to  the  simple  step- 

pulse  form  and  can  be  provided  in  any  given  case  by  a crash  analysis 

such  as  IITRAIN  or  for  parametric  studies,  the  values  of  V . . and 

01  j 

aoij  can  ke  vari-ec*  as  desired. 

Activity  4 represents  input  data  relevant  to  anthropometric 
parameter  values  such  as  mean  weights  of  a particular  body 
component  and  related  distribution  characteristics.  The  specific 
values  used  in  each  of  the  basic  passenger  configurations  will  be 
discussed  in  the  appropriate  section  of  this  chapter  as  they  are 
required.  In  examining  the  influence  that  the  response  character- 
istics of  the  surface/structure  at  the  impact  point  has  upon  the 
levels  of  injury  produced  it  became  clear  that  the  value  of  the 
impact  velocity  was  significant  and  that  care  should  be  taken  to 
establish  the  value  as  accurately  as  is  feasible  within  the  limits 
of  the  known  details  of  the  event.  Furthermore,  the  analysis  of 
passenger  motions  indicated  that  variances  in  the  impact  velocity 
would  occur  due  to  anthropometric  variances  and  other  uncertainty 
factors  such  as  slip  and  passenger  alertness..  For  these  reasons 
the  specific  passenger  population  in  any  basic  event  (that  is  the 
ith  consist,  jth  car,  and  kth  passenger  configuration)  was  further 
subdivided  into  10  equal  percentile  subpopulation  groups  and  the 
impact  velocity  assigned  to  each  subpopulation  group  was  given 
as 

v = v0  *i  ?l  (125) 

where 

VQ  = change  in  velocity  of  car 

^ = mean  value  of  dimensionless  impact  velocity 

« 

- distribution  factor  for  £th  subpopulation  group* 

In  this  manner  the  value  of  the  impact  velocity  was  allocated  in  a 
distributed  fashion  over  the  basic  population  to  appropriately 
reflect  the  totality  of  variances  which  will  exist.  The  anthropo- 
metric data,  in  many  instances,  appears  to  follow  a "normal" 
distribution  and  therefore,  the  "normal"  distribution  function  was 
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This 


assumed  to  apply  for  all  uncertainties  in  their  aggregate, 
function  is  illustrated  in  Figure  57  as  curve  A and  its  integrated 
or  accumulative  form  as  curve  B.  The  unit  of  standard  deviations 
was  used  in  this  representation  where  the  mean  (most  probable) 
value  corresponds  to  the  value  zero.  The  corresponding  incremented 
form  for  curve  B is  shown  (as  a dotted  line)  to  illustrate  the 
resolution  achieved  by  this  number  of  subpopulation  groups  and 

therefore  to  justify  the  selection  of  10  subdivisions.  The  general 

* 

approach  used  to  establish  the  values  of  is  illustrated  in 

Figure  58.  The  standard  deviations  to  mean  value  ratio  (A/gj)  was 
estimated  for  the  dimensionless  impact  velocity  functions,  for 

both  the  forward  facing  and  sideward  facing  passenger  configura- 
tions. Since  the  ratios  for  each  of  these  configurations  were 

* 

similar  a single  set  of  values  were  selected  for  use  in  the 

methodology.  Furthermore,  since  the  factor  is  symmetric  about  the 
mean  value  only  five  specific  values  needed  to  be  specified.  These 
five  are  given  in  the  "fixed"  data  section  as  PF(1)  to  PF(5).  The 
remaining  values  (to  PF(10))  are  computed  in  the  preliminary  calcu- 
lation section  of  the  code. 

Activity  5 represents  the  introduction  of  passenger  spacing 
details  and  other  means  for  dealing  with  anthropometric  variations 
or  the  likelihood  of  reducing  or  eliminating  the  injury.  These 
details  will  be  presented  in  the  corresponding  section  of  this 
chapter  for  each  of  the  basic  passenger  configurations. 

Activity  6 corresponds  to  the  establishment  of  the  impact 
velocity  or  eQuivalent  applicable  to  the  number  of 

passengers  Nj_jk£(m)  in  the  £th  subpopulation  group  involving  the 
mth  body  component.  Different  procedures  were  used  in  each  of  the 
basic  passenger  configurations. 

Activity  7 consists  of  the  introduction  of  threshold  values 
SIC  separating  each  of  the  adjacent  injury  level  categories  for  the 
severity  index.  The  values  for  the  head,  face  and  chest  as  indi- 
cated in  Figure  3 were  used  and  are  included  in  the  "fixed"  data 
section  of  the  code.  Whiplash  type  neck  injuries  were  arbitrarily 
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FIGURE  57.  NORMAL  DISTRIBUTION  FUNCTION 


o 50  100 

Population  Percentile  Group 

FIGURE  58.  DISTRIBUTION  FACTOR  FOR  POPULATION  GROUPS 
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limited  to  the  two  lower  injury  levels  of  minor  and  moderate.  The 
critical  acceleration  (see  equation  (40))  of  12  g was  selected  as 
the  threshold  value  separating  "no  injury”  from  "minor  injury”  and 
a corresponding  value  of  20  g was  selected  for  the  minor/moderate 
boundary.  An  upper  limit  control  value  of  100  was  used  to  limit  the 
injury  levels  to  "moderate".  These  data  are  also  included  in  the 
fixed  data  section  of  the  code. 

Activity  8 consists  of  the.  introduction  of  data  for  the 
properties  of  the  interior  components  of  the  car.  It  has  been 

assumed  that  each  car  in  each  consist  is  identical  thus  this  data 
group  is  "input"  for  each  consist.  Four  types  of  interior 

components  are  treated  and  defined  as  follows  by  a component  index 
JR  where 

JR  = 1 defines  a panel 
JR  = 2 defines  a seatback 
JR  = 3 defines  a stanchion 
JR  = 4 defines  a floor 


Furthermore,  each  component  type  is  assigned  a response  index,  IR, 
where 

IR  = 1 defines  a perfect  absorber 
IR  = 2 defines  an  elastic  response 
IR  = 3 defines  an  inelastic  response 

Two  response  constants,  RC1  and  RC2,  are  needed  to  define  the 
response  of  any  given  impact  situation.  RC2  is  only  required  for 
the  inelastic  response  case.  These  data  are  user  input  on  a problem 
to  problem  basis. 


Activity  9 is  that  activity  which  computes  the  value  of  the 
severity  index  corresponding  to  the  sub-subpopulation  group  for 
which  the  stimuli  have  been  determined.  A subroutine  IMPACT 
performs  this  task  given  the  impact  value,  the  structure  type,  and 
the  corresponding  structure  parameter  values.  Two  anthropometric 
constants  AH  and  WH  corresponding  to  the  contact  area  and  the 
weight  of  the  head  are  used  in  this  subroutine.  These  constants  are 
introduced  into  the  fixed  data  section  of  the  code. 
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Activity  10  is  that  activity  which  constructs  the  passenger 
injury  profiles.  A subroutine  COMPAR  compares  the  computed 
severity  index  SIA  for  each  sub-subpopulation  group  with  the 
critical  values,  SIC,  for  the  corresponding  mth  body  component  and 
accumulates  the  number  of  passengers  which  sustain  a given  level  of 
injury.  Since  the  population  group  is  generally  small  the  number 
of  passengers  which  appear  in  any  injury  level  may  appear  as  a 
decimal  quantity.  The  proper  interpretation  of  this  number  is  the 
probabilistic  interpretation.  Thus  when  the  results  indicate  that 
10.8  passengers  (out  of  a total  of  631)  sustain  an  injury  of 
"moderate  intensity"  it  means  that  the  probability  of  sustaining 
such  an  injury  is  0.017  for  the  totality  of  details  prescribed  by 
the  problem.  The  passenger  injury  profiles  are  represented  in 
terms  of  the  number  of  passengers  sustaining  a given  level  of 
injury  for  each  kth  passenger  configuration,  in  each  jth  car  of 
each  ith  consist.  These  values  are  totaled  twice  to  yield  a break- 
down by  car  and  a summary  for  the  entire  crash.  While  it  may  be  of 
interest  to  collect  and  print  out  more  detailed  injury  profiles, 
such  as  by  body  component,  or  to  express  the  profiles  in  terms  of 
percents  or  probabilities  or  other  entities,  such  an  output  format 
was  not  implemented  at  this  time. 

Activity  11  represents  an  evaluation  of  the  structural  crash- 
worthiness  of  the  cars  involved  in  the  crash  and  is  outside  the 
scope  of  this  passenger  injury  prediction  methodology  task. 

5.2  Standing  Passenger 

The  motion  details  of  the  standing  passenger  were  developed 
and  evaluated  in  Section  4.4.  A commonalty  of  behavior  was  found 
which  indicated  that  if  an  impact  with  the  floor  occurred  as  a 
result  of  the  motion  that  the  impact  velocity  would  be  independent 
for  all  conditions  examined  (both  car  motion  and  anthropometric 
parameters).  Thus  a critical  floor  impact  velocity,  VIFC  of 
22.4  ft;  was  selected  and  introduced  in  the  fixed  data  section  of 
the  code  for  this  effect.  The  commonalty  of  behavior  also 
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indicated  that  when  impact  with  a barrier  (panel  or  stanchion) 
occurred  that  the  impact  velocity  for  those  parameters  examined  is 
equal  to  the  sum  of  the  initial  car  velocity  VQ  and  a velocity 
increment  AVX»  The  latter  velocity  increment  depends  upon  the 
"time  of  flight"  of  the  passenger;  the  time  of  flight  was 
approximated  by  dividing  the  distance  xb  between  the  passenger  and 
the  barrier  by  the  car  velocity  V . This  ratio,  now  expressed  as 

ZT  = ^ = Z/VO  (126) 


is  used  as  the  independent  variable.  The  results  shown  in 
Figure  26  were  expressed  in  an  analytical  form  over  two  ranges  of 
ZT  as  shown  in  Figure  59.  These  approximations  are 


X 

> 

<3 

= B1 

+ D 1 (ZT)  , 

0 < ZT  <_  Cl 

and 

iVx 

= A 1 

- G 1 (ZT  - 

Cl  )2. 

Cl  < ZT  < El 

( 127) 

where 

B1  = 

A1 

- Cl  ( D 1 ) 

and 

G 1 = 

(FI 

- A 1 ) / ( El  - 

Cl)2. 

The 

constants 

A1  , Cl,  D 1 , 

El , and 

FI  are  incorporated  into  the  fixed 

data 

section  of  the  code. 

If  ZT  is  greater  than  El  then 

impact  with 

the  floor  will  occur  at  the  impact  velocity  VIFC. 

Figure  60  presents  the  program  details  for  the  standing  pas- 
senger. The  value  of  the  index  K is  1 for  this  passenger  configura- 
tion and  these  steps  are  taken  for  the  population  group  NPG(I,J,1), 
that  is  for  each  ith  consist  and  each  jth  car.  XF  is  the  total 
number  of  passengers  impacting  upon  the  floor,'  hence,  statement  133 
represents  an  initializing  step.  The  population  is  divided  into  10 
subpopulation  groups  equal  to  XL  in  number.  These  passengers  will 
be  distributed  equally  spacewise  between  the  barrier  or  stanchion 
spacing  SF  (input  data)  for  each  ith  consist.  Statements  135  to  157 
are  repeated  for  each  of  these  subpopulation  groups.  The  spacing  Z 
and  the  associated  time  ZT  are  then  computed.  IF  ZT  is  greater 
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FIGURE  59.  VELOCITY  INCREMENT  FOR  STANDING  PASSENGER 
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FIGURE  60.  PROGRAM  DETAILS  FOR  STANDING  PASSENGER 


than  E1;  floor  impact  will  occur  and  the  population  group  is  added  to 
the  accumulator  XF.  If  ZT  is  less  than  El,  the  proper  barrier  impact 
velocity  VI  is  computed.  The  effect  of  impact  occurring  while  the 
barrier  is  still  in  motion  included  by  letting  the  VO  contribution 
be  the  minimum  between  VO  and  VO*SQRT( Z/XO ) where  XO  is  the  appro- 
priate car  stopping  distance;  this  follows  from  equation  (51). 

A distribution  factor  DF  is  then  applied  which  further 
subdivides  the  subpopulation  group  into  two  groups  reflecting  the 
fraction  DF  of  the  population  which  is  exposed  to  a panel  and  the 
fraction  (IDF)  which  is  exposed  to  a stanchion.  Fifty  percent  of 
those  exposed  to  a stanchion  are  assumed  to  impact  the  stanchion; 
the  other  half  are  assumed  to  miss  the  adjacent  stanchion  and  are 
assumed  to  impact  the  floor  regardless  of  the  closeness  of  the  next 
stanchion.  Finally,  to  account  in  some  manner  for  passenger  alert- 
ness and  the  fact  that  each  of  the  above  impacts  may  occur  at  body 
positions  which  are  relatively  invulnerable  to  the  impact  process, 
such  as  the  arms,  the  population  was  split  into  three  segments.  The 
factor  PH ( 1 ) suffered  head  impact,  the  factor  PH(2)  suffered  chest 
impact  and  the  remaining  are  uninjured. 

Statements  144  to  154  performed  the  appropriate  population 
accounting,  calls  IMPACT  to  compute  the  actual  severity  index  SIA 
and  then  calls  COMPAR  to  evaluate  the  corresponding  level  of  injury 
which  has  occurred  for  each  sub-subpopulation  group  which  impacts  a 
barrier.  The  code  then  accumulates  the  number  of  passengers 
sustaining  a given  level  of  injury  in  the  array  XIG ( I , J, K , I J) . 
Statements  158  to  164  then  repeat  this  process  for  floor  impacts. 
Floor  impacts  which  occur  at  low  speeds  (i.e.,  low  values  of  VO)  are 
reduced  in  magnitude  by  selecting  the  minimum  between  VO  and  VIFC. 

5.3  Forward  Facing  Seated  Passenger 

The  motion  details  of  the  forward  facing  seated  passenger  were 
developed  and  evaluated  in  Section  4.5.  That  analysis  showed  that 
the  impact  velocity  of  the  passenger  with  the  object  immediately 
forward  of  the  passenger  was  a function  of  stopping  distance  param- 
eter, k,  and  that  a range  of  uncertainties  occurred  as  a result  of 
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anthropometric  variances  and  the  degree  of  hip  loft  and  seat 
slippage  which  could  exist.  A mean  value  curve  ef.  was  constructed 
as  illustrated  in  Figure  61  and  which  consisted  of  two  analytical 
parts.  For  small  values  of  ic(the  symbol  C is  used  in  the  code)  a 
simple  square  root  dependence  is  used.  For  values  of  k greater  than 
C2}  a linear  relationship  involving  the  costants  B2  and  D2  is  used. 
The  intercept  value  B2  is  obtained  from  A2,  C2  and  D2  in  the  pre- 
liminary calculation  section  of  the  code.  A measure  of  the  vari- 
ance (assumed  to  be  equal  to  one  standard  deviation)  is  shown  in 
this  figure. 

Figure  62  presents  the  program  details  for  the  forward  facing 
seated  passenger.  The  value  of  the  index  K is  2 for  this  passenger 
configuration  and  these  steps  are  applied  to  the  population  group 
NPG(I,J,2).  The  stopping  distance  parameter,  C,  is  evaluated  using 
the  corresponding  body  dimension  BD  (for  the  c.g.  location)  and  the 
stopping  distance  XO.  The  stopping  distance  is  evaluated  at  the 
beginning  of  the  basic  calculation  section  of  the  code  from  corre- 
sponding velocity  change  and  average  acceleration  input  data.  The 
step-pulse  acceleration  pulse  shape  is  used.  The  dimensionless 
mean  value  impact  velocity  FE  is  then  evaluated.  Statements  169  to 
184  are  then  repeated  for  each  Lth  percentile  group.  The  impact 
velocity  is  computed  first  from  equation  (125).  The  subpopulation 
is  then  divided  into  two  categories  by  use  of  the  distribution 
factor  DF.  The  DF  faction  is  subjected  to  impact  upon  panels;  the 
remainder  are  subject  to  impacts  upon  seatbacks.  The  subroutine 
IMPACT  is  then  called  to  obtain  a reference  severity  index,  SI. 
Statements  176  to  184  are  then  repeated  for  each  of  the  three  body 
components  (head,  face,  and  chest).  If  a panel  (JR=1)  is  involved^ 
only  head  injuries  are  considered  and  subroutine  COMPAR  is  called 
to  evaluate  the  level  of  injury  involved.  The  final  step  (state- 
ment 183)  is  the  accumulation  of  the  number  of  passengers  involved 
in  the  resulting  injury  level  (IJ). 

♦ 

A somewhat  more  involved  procedure  is  used  if  a seatback 
(JR=2)  is  involved.  This  is  done  in  an  attempt  to  properly  dis- 
tribute the  subpopulation  over  the  body  components  which  could  be 
involved.  This  procedure  is  illustrated  in  Figure  63* 
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FIGURE  61.  MEAN  IMPACT  VELOCITY  FUNCTION  FOR  FORWARD 
FACING  SEATED  PASSENGER 
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FIGURE  62.  PROGRAM  DETAILS  FOR  FORWARD  FACING  SEATED  PASSENGER 
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FIGURE  63.  DISTRIBUTION  OF  BODY  COMPONENT  INVOLVED  IN  IMPACT 


A trilinear  torso  length  distribution  curve  for  the  population  was 
selected  and  is  characterized  by  the  two  constants  BH1  and  BH2.  A 
cord  length,  SC,  (evaluated  in  the  preliminary  calculation  section 
of  the  code)  was  then  computed  from  the  seatback  height  SH  (above 
the  hip)  and  the  space  SS  between  the  passenger  c.g.  and  the 
seatback  immediately  forward  of  the  passenger.  The  faction  of  the 
population  (1-PP(1))  whose  torso  length  is  less  than  the  cord 
length  will  be  exposed  to  a head  impact.  The  faction  of  the 
population  (PP(  1 )-PP(2) ) whose  torso  is  less  than  the  sum  of  the 
cord  length  and  the  face  length,  DS,  will  be  exposed  to  a face 
impact.  The  remainder  will  be  exposed  to  chest  impacts.  The  actual 
severity  index  SIA  is  computed  from  the  reference  severity  index  SI 
by  using  a factor  AF.  This  factor  was  intended  to  provide  some 
modification  of  the  severity  index  due  to  the  change  in  weight  of 
the  body  component  involved.  This  feature  was  not  incorporated 
into  the  current  methodology,  therefore,  a value  of  AF  of  unity  is 
used.  The  parameters  A2,  C2,  D2,  BD,  DS,  BH1  and  BH2  are  incorpo- 
rated into  the  fixed  data  section  of  the  code. 

5.4  Sideward  Facing  Seated  Passenger 

The  motion  details  of  the  sideward  facing  seated  passenger 
were  developed  and  evaluated  in  Section  4.7.  That  analysis  showed 
that  the  impact  velocity  of  the  passenger  with  the  object  immedi- 
ately to  his  or  her  side  was  a function  of  the  stopping  distance 
parameter,  k,  and  that  a range  of  uncertainty  occurred  as  a result 
of  anthropometric  variances  and  spacing  from  the  adjacent  object. 
A mean  value  curve  was  constructed  as  illustrated  in  Figure  64 
which  consisted  of  two  analytical  parts.  For  small  values  of  k 
(the  symbol  C is  used  in  the  code)  a simple  square  root  dependence 
is  used.  For  values  of  k greater  than  C3  a constant  value  A3  was 
used  for  the  value  of  gj.  . A measure  of  the  variance  (assumed  to  be 
equal  to  one  standard  deviation)  is  shown  in  this  figure. 

Figure  65  presents  the  program  details  for  the  sideward  facing 
seated  passenger.  The  value  of  the  index  K is  3 for  this  passenger 
configuration  and  these  steps  are  applied  to  the  population  group 
NPG( I , J, 3) • 
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FIGURE  64.  MEAN  IMPACT  VELOCITY  FUNCTION  FOR  SIDEWARD 
FACING  SEATED  PASSENGER 
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FIGURE  65.  PROGRAM  DETAILS  FOR  SIDEWARD  AND  BACKWARD  FACING  SEATED  PASSENGER 


The  stopping  distance  parameter,  C,  is  evaluated  using  the  corre- 
sponding body  dimension,  BE  (for  the  c.g.  location)  and  the 
stopping  distance  XO.  The  dimensionless  mean  value  impact  velocity 
FE  is  then  evaluated. 

The  number  of  passengers  who  will  suffer  an  impact  will  be 
limited  to  those  at  the  one  end  of  the  side  bench  (or  individual 
seats)  that  is  those  adjacent  to  an  obstacle.  Thus  a reduction 
factor  SP  is  used  which  is  equal  to  one-half  (for  one  end  of  the 
bench)  divided  by  the  number  of  passengers  on  the  side  bench.  This 
number  is  obtained  by  dividing  the  bench  length,  SB,  by  the  body 
width,  H3«  The  remaining  procedures  are  similar  to  those  used  for 
the  forward  facing  seated  passenger.  The  values  of  A3,  C3,  and  H3 
are  incorporated  into  the  code  in  the  fixed  data  section  of  the 
code . 

5.  5 Backward  Facing  Seated  Passenger 

Passengers  in  this  passenger  configuration  will  only  be 
subjected  to  a "whiplash”  type  of  injury  which  is  largely  dependent 
upon  the  level  of  acceleration  (see  Section  4.6).  The  program 
details  for  this  passenger  configuration  are  presented  in  Fig- 
ure 65.  The  index  K=4  for  this  case  and  the  computational  steps  are 
applied  to  the  population  group  NPG(I,J,4).  The  population  is 
divided  into  10  percentile  subgroups  according  to  the  index  L and 
the  car  acceleration  AO  is  adjusted  by  the  factor  PF  in  each  Lth 
case  to  obtain  some  variance  due  to  sensitivity  and  anthropometric 
variations.  These  accelerations  are  then  compared  to  the  criteria 
values  with  the  aid  of  subroutine  COMPAR  and  the  accounting  proce- 
dure with  respect  to  the  number  of  passengers  sustaining  a given 
level  of  injury  is  followed  as  was  done  with  the  other  passenger 
configurations . 

5.6  Typical  Results 

The  results  of  a hypothetical  crash  situation  are  presented  to 
demonstrate  the  nature  of  the  output  of  the  methodology.  The  crash 
environments,  i.e.,  the  change  in  velocity  VO  and  average  accelera- 
tion AO  for  each  car  of  each  consist,  were  selected  arbitrarily  as 
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were  the  car  interior  component  properties,  spacing  and  other 
factors.  Hopefully,  these  values  are  somewhat  realistic  such  that 
the  numerical  results  will  be  meaningful.  The  car  interior 
arrangements  and  the  properties  of  some  of  the  interior  components 
were  selected  to  demonstrate  some  of  the  influence  that  these 
factors  have  upon  the  population  of  injuries  in  an  intracity  train 
crash  situation. 

The  results  are  included  at  the  end  of  this  section.  A summary 
is  presented  of  the  number  of  passengers  involved  in  the  crash. 
This  is  followed  by  a detailed  breakdown  of  the  passengers  by 
number  with  respect  to  passenger  configuration,  car  and  consist. 
In  this  example,  consist  number  1 did  not  have  any  sideward  facing 
seats.  The  cars  of  the  second  consist  are  more  conventional  in  this 
respect.  The  first  consist  is  rather  crowded  and  therefore,  there 
is  a larger  number  of  standing  passengers  than  in  the  second 
consist. 

The  injury  analysis  summary  comprises  the  total  number  of 
passengers  in  each  of  the  six  injury  categories  and  a corresponding 
breakdown  by  car  and  consist.  The  corresponding  details  with 
respect  to  each  passenger  configuration  are  also  presented.  The 
remaining  output,  including  the  "fixed  data"  for  this  example,  are 
presented  in  the  appendix. 

Most  of  the  fatalities  which  occur  in  this  example  are  the 
result  of  injuries  sustained  by  standing  passengers  in  consist 
number  1 who  apparently  fall  to  the  floor.  The  floor  system  for  the 
cars  in  this  consist  is  quite  stiff  whereas  the  stanchions  and 
panels  are  of  nominal  stiffness.  The  distribution  of  the  injuries 
sustained  by  standing  passengers  in  consist  number  2 is  quite 
different  and  reflects  somewhat  stiffer  stanchions  and  panels  and  a 
considerably  softer  floor. 

The  converse  is  true  for  the  forward  facing  seated  passenger. 
The  seatback  properties  in  the  cars  of  consist  number  2 are  stiffer 
and  result  in  a definite  shift  in  the  intensity  of  the  injuries  in 
this  configuration.  Only  a few  minor  injuries  are  sustained  by 
passengers  in  this  configuration  for  consist  number  1 even  though 
many  more  passengers  are  involved. 
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No  comparison  can  be  made  for  the  sideward  seated  passengers 
since  none  existed  in  consist  number  1.  The  number  of  injuries 
sustained  by  sideward  seated  passengers  in  consist  number  2 is 
quite  small  and  is  due  in  part  to  the  fact  that  relatively  few  of 
these  passengers  are  close  to  an  interior  object.  Finally,  no 
whiplash  injuries  occur  in  the  backward  facing  seated  passenger 
configuration  in  either  consist  due  to  the  relatively  low  accelera- 
tion levels  involved. 


Intra-city  train  crash  injury  analysis 
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INJUHY  ANALYSIS  summawy 
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APPENDIX 


COMPUTER  CODE  DESCRIPTION 


A brief  descr  iption  of  the  computer  code  used  in  the  injury 
prediction  methodology  entitled  "Intra-City  Crash  Injury  Analysis" 
is  presented.  This  code  was  written  in  FORTRAN  V and  has  been  used 
on  a UNIVAC  1108  computer.  It  is  a relatively  small  code  and  has 
minimal  storage  requirements.  The  code  is  divided  into  eight 
sections . 

1 . Data  Section 

The  multitude  of  parameter  values  which  are  fixed  in 
value  are  presented.  The  parameters  include  dis- 
tribution values,  anthropometric  data,  injury 
criteria,  etc.  Note  that  as  some  of  these  data 
values  are  determined  more  precisely,  this  section 
may  be  modified. 

2 . Input  Section 

The  introduction  of  the  necessary  input  data  is 
provided.  These  data  and  their  corresponding  format 
requirements  are  presented  in  Table  A-1.  The  code 
has  a multiple  problem  solving  capability  such  that 
any  number  of  problems  can  be  run  by  adding  data 
card  groups  as  specified  in  Table  A-1. 

3 . Initialize  and  Preliminary  Calculations 

Certain  arrays  and  variables  are  initialized  and  a 
series  of  preliminary  calculations  which  facilitate 
the  subsequent  calculations  are  performed. 

4 . Printout  - First  Phase 

This  phase  prints  a title  and  problem  identification 
along  with  the  precrash  speed,  a breakdown  of  the 
total  number  of  passengers  involved  in  the  crash  and 
the  total  number  of  passengers  involved  in  each 
consist.  If  only  one  consist  is  involved,  this  fact 
is  printed.  The  passenger  number  distribution  by 
consist  and  by  car  for  each  of  the  four  passenger 
configurations  treated  in  the  methodology,  and  the 
corresponding  totals  and  subtotals  are  presented. 

5 . Basic  Calculations 

All  of  the  basic  calculations  are  performed  for  each 
of  the  four  passenger  configurations  as  described  in 
Section  5 of  this  report.  These  calculations  are 
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TABLE  A-l  INPUT  FORMAT 


Card  (Format) 

Variables 

Description  (units) 

Card  1 

TITLE 

Title 

(18A4) 

Card  2 

SPEED 

Relative  speed  of  consists  at  impact 

(F10. 1,4110) 

IC 

Number  of  consists  (<2) 

NC(1) 

Number  of  cars  in  Consist  1 (s<10) 

NC  (2) 

Number  of  cars  in  Consist  2 (v<10) 

IPP 

Print  control  index  for  stored  data 

Card  Group  A (one  group  for  each  consist  1=1, IC) 

Card  Group  A1  (one  subgroup  for  each  kth  passenger  configuration,  K=l,4) 


Card  Al-1 
(1015) 

Card  Al-2 
(10F5.2) 


K=1  Standing  passenger 
K=2  Forward  facing  seated  passenger 
K=3  Sideward  facing  seated  passenger 
K=4  Backward  facing  seated  passenger 


NPG(I,J,K)  Number  of  passengers  in  kth  configuration 
J=1,NC(I)  of  jth  car  of  ith  consist 


DF(I,J,K)  Distribution  factor  for  kth  passenger 

J=1,NC(I)  configuration  of  jth  car  of  ith  consist 


NOTE:  The  last  card  of  this  group  corresponding  to  Card  Al-2  for  K=4  is  blank 

since  DF(I,J,4)  is  not  used. 

(End  of  Card  Group  Al) 


Card  A2 
(10F8.2) 

Card  A3 
(10F8.2) 

Card  A4 
(4110) 


Card  A5 
(4E10.4) 


Card  A6 
(4E10.4) 


Card  A7 
(8F10.2) 


AO ( I , J ) 
J=1,NC(I) 

Acceleration  of  jth  car  of  ith  consist (g) 

V0(I,J) 

J=1,NC(I) 

Velocity  change  of  jth  car  of  ith  consist 
(fps) 

IR (I , JR) 
JR=1 ,4 

Response  Form  Index 
IR=1  for  perfect  absorber 
IR=2  for  elastic  response 
IR=3  for  inelastic  response 

RC1 (I , JR) 
JR=1 , 4 

Response  constant  Number  1 defined  as 

Yield  stress  (psi)  for  IR=1 

Spring  constant  (lb/in.)  for  IR=2  and  3 

RC2 (I , JR) 
JR=1 , 4 

Response  constant  Number  2 defined  as 
Critical  force  (lb)  for  IR=3  otherwise 
not  used  (set=0.) 

SH(I) 

Height  of  seatback  above  hip  (ft) 

SS(I) 

Space  between  forward  facing  seated  pas- 
senger and  object  in  front  of  seat  (ft) 

SB(I) 

Side  bench  length  (ft) 

SF(I) 

Floor  space  length  (or  stanchion  spacing) (ft) 

For  additional  problems  repeat  entire  group  (Cards  1 to  A7). 
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for  each  car  in  each  consist.  Two  subroutines  are 
called  from  this  section.  The  first,  subroutine 
IMPACT;  computes  the  severity  index  in  accordance 
with  the  appropriate  surface/structure  response 
form.  The  second,  subroutine  COMPAR,  compares  the 
value  of  the  severity  index  with  the  appropriate 
critical  values  and  determines  the  level  of  injury 
sustained.  The  number  of  passengers  associated  with 
this  injury  level  are  accumulated  for  each  passenger 
configuration,  for  each  car,  and  for  each  consist 
involved  in  the  crash  event. 

6 . Printout  - Second  Phase 

This  phase  prints  out  the  injury  data.  A summary  is 
presented  first  which  lists  the  total  number  of 
passengers  assigned  to  each  injury  level  including 
uninjured  together  with  a corresponding  breakdown  by 
car  and  consist  including  appropriate  subtotals. 
This  is  followed  by  injury  profiles  for  each  of  the 
passenger  configurations.  These  details  are  given 
for  each  car  in  each  consist. 

7.  Printout  - Third  Phase 

This  phase  prints  out  the  remaining  input  data  to 
document  the  problem  appropriately.  This  includes 
the  interior  components  data  for  each  consist,  that 
is  the  response  from  index  (IR)  and  the  two  associ- 
ated response  constants  RC1  and  rC2  for  each  of  the 
four  interior  component  categories  (JR).  It 
includes  the  car  crash  environments  (velocity  change 
and  acceleration)  for  each  car  of  each  consist  and  a 
collection  of  distribution  factors  and  spacing 
parameter  values.  The  former  DF  is  presented  for 
three  of  the  four  passenger  consists  and  for  each 
car  of  each  consist.  The  latter  (SH,  SS,  SB,  and 
SF)  are  presented  for  each  consist. 

8 . Printout  - Fourth  Phase 

This  phase  is  an  optional  output  group  controlled  by 
the  control  index  IPP  (IPP  >0  provides  for  the 
printing  of  this  information).  These  data  are  the 
so-called  ’’fixed  data”  used  in  methodology  and 
include  probability  distribution  factors  PF(I) 
(1=1,5)  and  a collection  of  constants  used  to  define 
the  results  of  the  motion  analyses.  These  data  also 
include  some  anthropometric  data  and  the  severity 
index  criteria  data.  The  optional  form  is  intended 
to  provide  periodic  documentation  of  the  fixed 
values  under  conditions  where  some  of  these  values 
may  be  updated  or  modified. 
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A detailed  listing  of  the  code  follows  together  with  a listing 
of  the  two  subroutines  IMPACT  and  COMPAR.  A sample  problem  was  runj 
and  the  corresponding  output  data  are  presented  in  Section  5 of  the 
report  with  the  exception  of  the  output  from  Sections  7 and  8 of 
the  code  which  are  presented  on  the  following  pages  of  this 
appendix  and  correspond  to  the  sample  problem. 
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IMKA-C1TY  TRAIN  CRASH  I N J U K Y ANALYSIS 

DIMENSION  inLE(l*)*AG(2*10)*VU(2*10)*AO(2*10)*NCl2)»NPTl2)* 
1NPC(2*10)*NPG(?*10*P)*NPGTC(2*P)*NPGlHP)*PF(10)*XITuTl6J*SP(2)* 
2xIT(2.6)*XlC(2*l0*e)*S6(2)*DF(2*10*P)*6H(2)*SS(2)*SC(2)*PPl2)» 
3Pm(2*3)*XIG(2*1G*p*6)*PH(2)*SF(2) 

CORMON/SUpKl/  1*J*R*VI*AH*RH,1R(2*P)  ♦ R C 1 ( 2 ♦ u ) »RC2(2*P) 

1 /SICKIT/  SK(P*S) 


DAI  a 5ECTIGN 


DATA 

DATA 

1 22. a 

DATa 

1 

DATA 

2 

DATA 

3 

DAT  A 
DATA 
DATA 


PF  ( 1 ) 
V IFC  » 
3.6*. 

S I 0 ( i 
6 5 0 • ♦ 
SIC  (2 
buO  • * 
S I C ( 3 
1 0 0 0 . 
S1CIP 

p D * A h 

Pn(  1 ) 


♦ PF ( 2 ) ♦PFC3) 
A 1 • C 1 • 0 1 * E 1 * 
35*1  1 .03* .71 

* 1 ) *SIC(1 *2) 
1000. *1333. f 
» 1 ) *SIC ( 2*2) 
1000,* 1500** 

♦ 1 ) ♦ S T C ( 3 ♦ 2 1 
*1500.*200u, 
*l)*SlC(P*2) 
*WM*t?Hl*rtH2* 
*Ph(2) /,P0* . 


,PF(P).PF(b)/.72*.8?*.66*.9P*.98/ 

F 1 . A3*C3 * A2.C2*  D2*m3/ 

» -P • * 1 • 3 1 * , Obb  * 1 ,7u*,bb*,1 *1,75/ 
•SICC1  *3)  *S1C11*P)  *S1C(1*5)/ 

1 6 b 7 . * 2 0 0 0 . / 

*SiC(?*3)*SIC(2*a)*SlC(2*S)/ 

2f  00  , * 2500  ,/ 

*SIC(3*3) *SIC(3*p) * S I C ( 3 ♦ b ) / 

*2500. *30  0 0./ 
*S1C(P*3)/1?.*20.*100./ 
DS*6t/l.*10.*10.*1.50*2.50*,7bf.27/ 
10/ 


C 

C 

C 


Input  section 


10  -READ  (5* 100* END  = 9P9)  T 1 T L E * SP E E D * I C * nC C 1 ) . NC ( 2 ) * I PP 
luO  FORMAT  ( 1 8A4/F 1 0 . 1 * a 1 1 0 ) 

DU  IIP  Is  1 * I C 

- JC  = NC(I) 

U 0 1 0 2 K s 1 « U 

READ  (b*10u)  (NPG( I* J.K) * Jsl * JC) 

READ  (5*105)  U)F(I*J*K) *J=1*JC) 

F Ohm  AT  (1015) 

FORMAT  (10F5.2) 
hFAO  (5*108)  ( A0( I * J) . J=1 * JC) 

READ  (5*108)  CVGU*J)*J=1*JC) 

FORMAT  (10F8.2) 

READ  ( 5 ♦ 1 1 u ) UR(I*JR)*Jksi,<j) 

110  F Ohm  AT  («I10) 

READ  (5*112)  (RC  1 ( I * JR) ♦ JRsl ,p) 

112  Format  (PF10.P) 

PF AO  (5*112)  ( RC2 ( I • JR ) * JR= 1 , P ) 

114  READ  (5*116)  Sh( I ) *SS( I) *Sb( 1 ) *SF ( I) 

116  FORMAT  (8F10.2) 


1 u2 
1 U P 

lo  5 

1 0 6 
lob 
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C 

0 INIUAUZP  A,\r  HP£L  IMNAKY  CALCULATION 
L 

N P T 0 T = U 

L'Q  POO  1 = 1 t 1C 

JC  = N'C(l) 

N P 1 ( 1 ) = 0 
5 P ( I ) = .S*H3/SHtI) 

been  = s u k t c s h ( i ) * s h c i ) f s s ( i ) * s o ( n ) 

DO  P O o J = If J C 

p o o n p c i n j ) = o 


2 uP 
?.u  4 
p U h 


?0  8 

PI  0 

PI? 


214 


21b 


21  K 


UP  POh  1 = 1 t IC 

JC  = ^C(I) 
i)u  20  4 J=1»JC 

0 U 2 u P * = 1 t 4 

NPCIItJJ  = NPC(I*J)+NPG(ItJfK) 

N p T C I ) = NPT  UH'JPCC  I t J) 
uPTuT  = iviPTuT  + NPl  ( n 
GO  210  K:ha 
NPb T T ( * ) = 0 

00  210  1=  1 tlC 

JC  = f'  C ( I) 

N P G T L ( 1 t K ) = 0 
on  20b  J=lfJC 

OPGTC  l l • * ) = ^PGTCCI  tK)+i\PG(ItJti\) 
N P b IT  ( *)  = NPGlTl^J+NPGTCdfK) 

0 0 21?  L = o t 1 0 

L L = 11  -L 
PF(L)  = 2.-PFCLU 
b?  = A2-C2*D2 


t*  1 = A 1 - C 1 * L>  1 


G 1 = ( F 1 - M ) / ( (El-Cl  ) **?) 
bL  = hbP-Bbil 
on  2 1 b I J=  i to 
X110HIJ)  = 0 . 


0 0 21b  I = 1 1 1 C 
JC  = N-C(l) 

X I 1 ( I t I J ) = o . 

on  2i4  jsitjc 

X I C ( 1 t J t 1 J ) = 0. 

00  214  KSlffi 
XIG(  1 • JtK, JJ)  = u • 

00  21b  b = 1 t 2 

PP(M)  = 1 .-SL»(SC( D+FLOAT (P-1 )*0S-Ph1 ) 
IF  (PP(«J  .LT*  0.)  PP(fi)  = 0. 

IF  (PP(b)  .GT . 1.)  PP(m)  = 1. 


CONTINUE 
PMlItl)  = 
P m ( 1 t 2 ) = 
H * ( I t 3 ) = 


1 .-PP ( 1 ) 

PP( 1 )-PP(2) 
P P ( 2 ) 
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L 

C PRIM  (JUT  - FIRST  PHASE  ( MJK  BEK  DISTkIbUTIUn  anP  I DF  n T 1 F L C A 1 i ON  ) 

L 

n P 1 T E (6 *ft00)  T I T LE * SPEFO* NPTOT 

boo  F OP'-'  M (lHlf2PX»U'MTPA-CITY  TP  A IN  CRASH  InJUKY  analysis'/// 
nox*  'SPECIFIC  C u m m E N T S - ' * 1 P a 4 / 

P 1 0 X * • k £ !_  A T J V F SPEFO  AT  CPASH  slfF6.1»'  «Hh'// 

32 ox* •total  N u m P E w OF  PA S SENG EPS  INVOLVED  In  C R a S h = ' * 1 3 ) 

IF  ( I C . t ft . 1 ) GO  TO  150 

ft  R I T E C ft  * 6 0 2 ) nPT ( 1 ) ,NP1 12) 

ft 02  FORMAT  (20x*» TOTAL  NUMBER  OF  PASSFNGEWS  IN  CONSIST  NO.  1 ='*16/ 
120X*«TuTAL  n U M B E ft  OF  PASSENGERS  IN  CuNSlbT  NO*  2 =»*lb///) 

GO  To  152 
ISO  * s 1 T F.  ( ft  » <-  0 u ) 

noa  FOft^AT  ( 2 ' X ♦ 1 l 0 n L Y ONE  CONSIST  INVOLVED)'///) 

1 52  * ft  I T E ( ft  * ft  1 0 ) 

ft  1 0 FORMAT  (///*3ox*  'PASSENGER  DISTPIBoTIOn  ft  Y n u R ft  F k • / / 

1 3 X » •CONSIST  CAP  total  NUMBER  S I A N I)  I N G SEATED  FACINo 

? SEATED  facing  SEATED  FACING' /sx* »Nu. • *7X* »nO, • *BX* I Iu  CAp'*22X 

3*  ' F O K -v  A P o • , B x » •SIDEWARD*  * H X * ' BACKWARD  '/) 

DO  1 b ft  1 = 1 tIC 
JC  = NC ( I ) 

L)U  1 bo  J=ltJC 

1 3 ft  R I T E ( ft  « ft  1 2 ) I*J*NPCCI*J)f(NPG(I*J*K)*Ksl*P) 

ft  1 2 format  (17*110*113*115*113*2116) 

1 5 ft  ft  p 1 T E C ft  * ft  1 p ) I*(nPGTC(I*k)*k=1*u) 

ft  1 a FQkmat  (/lx*  'TOTAL  In  CONSIST  NO  . ' ♦ I 2 * I 22 ♦ 1 1 3 * 2 I 1 ft// ) 

_ — _ iftKlTE  (fttftlft)  (NPGTT(k)*n=1»u) 

o 1 to  FORMAT  (//IX,  'TOTAL  IN  CftASH'  *130,113*21 lft) 
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C hASIC  CALCULATIONS 

c 

00  590  1 = 1 »IC 

JC  = NCC1) 
do  5pp  jsitjc 

xO(ltJ)  = ,5*V0( I tj) *V0( I »J) /( AU( I tJ) *32,2) 

K = 1 
XF  = 0. 

XL  = • 1 * F L U A T (ImPC»(I*J*K)  ) 

00  52  0 L = 1 t 1 0 
L = . 1 *SF ( I ) *FLOAT (L) 

ZT  = Z/VU(1»J) 

IF  (ZT  .GT,  El)  GO  TO  51b 

IF  (ZT  .GT.  Cl ) GO  Tu  510 

VI  = Vu(  I » J)*AMIN1  ( 1 .tSf>KT(Z/XO(  If  J)  ) )^bl+01*Zr 
GO  1 U 5 1 2 

510  VI  = V0( I * J ) * AMlul ( 1 . f SOkT ( z/ XU( 1 t J ) ) ) + A 1 + 01 ♦ ( ( Z T-C 1) **2) 
512  CONTINUE 

vO  510  JR=lf5*2 
X K = X L * 0 f C 1 t J f H ) 

IF  (JR  .GT.  1)  x P = .5* XL* ( 1 ,-UF  ( I » J » k ) ) 

IF  ( J K • G T • 1)  XF  = XF  + Xk 

XI&(I.J*Ktl)  = XIG(I*J*Kfl)+Xh»(l.-Ph(l)-PM(2)) 

CALL  IMPACT  ( J k . S I A ) 
no  510  MS  1,3* 2 
X M = Ph(  1 )*xk 

IF  (M  . C,  T . 1 ) X V s P H ( 2 ) * X K 

CALL  CUMPap ( si  At m.  IJ) 

510  XlG(ItJtKtlJ)  5 XlG(ItJ*KtIJ)+XM 
GO  Tu  520 
51b  XF  = XF+XL 
520  C 0 N T j n u E 

XIG(IfJfKtl)  = X 1 G C I t J t K , 1 ) + XF* ( 1 ,-PH( 1 )-Pn( 2) ) 

VI  = A M ] N 1 ( V I F c • v 0 ( 1 t J ) ) 

Call  1 mPac  T l o t S I A ) 

call  cg^papcsi At i , i j) 

XIG( 1 t JtK.  IJ)  s X IGC I t J . K t I J) +XF*Ph(  1 ) 

CALL  C0MPARlSlAt3t IJ) 

X IGC I t JtK t 1 J)  = Xl6(ItJtKtIJ)+XF*PMl2) 

K = 2 

C = .5*PD/XU(  I • J) 

FE  = A 2 * S Q k T ( C / C 2 ) 

I F (FE  .GT.  A2)  FE  = F2+D?*C 

UO  500  L=  1 t 1 0 

VI  = VQ(  ltJ)*F£*RF  (L) 
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DO  540  JK=1*2 

F 0 = 0 F ( I t J • K ) 

IF  (JR  .60.  2 ) FD  = l.-DF(I.J.K) 

XK  = . 1*FD*FL0AT (KPG( I » J*K) ) 

CALL  IMPACT  (JR. SI) 

DO  540  *=1.3 

IF  C Jr?  .tu.  1 .And.  m .GT.  1)  00  TU  54  u 

A F = 1 . 

SIA  = A F * S 1 

X M = P M ( I.M)*XR 

]F  (JR  . t «v  • 1)  X M = X R 

CALL  C 0 * P A * l 6 I A ♦ M f I J ) 

X 1 1,  ( I • J . K . I J ) = XJU(I.J.K.IJ)+XM 
54  0 L 0 N I I K' u F 
* = 3 

L = ,5*bt/X0(  1 » J) 

FE  = A3*SUP1 (C/Ci) 

IF  (FF  .GT.  A3)  FE  = Ai 

XL  = .l*FLUAT(NPG(I*JtK))*SP(I) 

x IG(  1 t J t K ♦ n = x 1 1-  ( T ♦ J ,*  t 1 ) + l 1 . - 5P  ( I ) ) ♦FLOAT  ( NPf,ll  ♦ J . K ) ) 

DU  55  0 L = 1 *30 

vl  = vu( 1 » J)*Ft*PF(L) 

00  550  JR=  1 » 2 
FD  = DF(I*J*K) 

IF  (JR  . t'  U . 2)  FU  = l.-DF(I»J»iO 

X k = X L ♦ F D 

CALL  IMPACT  (JR. SI) 

DO  550  Ms  1.3 

IF  (JP  .tO,  1 .and.  M ,GT.  1)  GO  TO  550 
AF  = 1. 

SIA  = A F * S I 
XM  = P M ( I . M ) ♦ X R 

1 F ( J P . 1 0 . 1 ) X M = X P 

CALL  CUMPAh  (SIA.M.IJ) 

XlG(IfJ»K*lJ)  S XlG(I.J»to.IJ)+X* 

550  CONTINUE 
k = a 

DO  570  Lsl  .10 
SIA  = A0(1»J)*PF(L) 

XM  = .1*FluAT(nPG(I.J.K)) 

CALL  C U m P a R ( S I a , a ♦ I J ) 

570  XIG(1 .J.k.IJ)  = XIG(1. J»K.IJ)+XM 
DO  5h4  I J si *6 
DU  5b 4 K s 1 . U 

564  XIC(I.J.IJ)  = XIC(I.J.IJ)+XIG(I.J.K.IJ) 

DO  5b A I JS 1 ♦ to  --  — 

5b  A XIKI.IJ)  = xn(I.IJ)+XIC(I.J.lJ) 

OU  590  I J=  1 . to 

590  XITOT(IJ)  = XITOT ( IJJ+XII  ( I.  1J) 
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C PRINT  OUT  - sec  UNO  Pm  ASF  (lNjiiKY  DATA) 

L 

vs  Wilt  ( vm  H ()  U ) (XiTlJlllJ1*IJsl»6) 

800  F 0 p '*  A T ( 1 H 1 ♦ <4  *5  X t < 1 «v  J U p Y ANALYSIS  SU^MkY'/// 

14bX.  'TOTAL  Un- INJURE U = 'tFe>.l/46X.  'TOTAL  Ml  NOP  s'  *Fo,l/ 

] 4bX  » I UT  *»L  w OOP  P AT  £ = ' tFo.l/abX.  ' T U T a L SERIOuS  s'  ♦ F (?  , 1 / 

2 4 6 X » 'TOTAL  CPITJCAL  = ' ♦ F b . 1 / 4 t X ♦ ' T u T A L FATALITIES  = »tFe>.l//) 
v*  k 1 T E l 8 » 6 u 2 ) 

8 o2  r 0 K P A T ( 4 1 X t ' 1 nJ.JP  Y PROFILES  8Y  LAP  A VO  CONSIST'// 

1 i X t ' L 0 N S 1 5 1 C A P ' t a 1 X t I I N J U K Y L t v F.  L ' / 5 X t ' N Q . • t 7 A t ' N U , ' . 7 X * ' u N 

2 N J 0 R t 0 ' ♦ 7 X » 1 i*  1 n 0 P ' » 7 x » 'MODERATE1  ♦ 7 x » 'SERIOUS'  *eX  t ' C H 1 T lCAL'  » 5 X * ' 

3 T * LI  T I & S ' / ) 

UO  810  1 = 1 tlC 

jc  = Nccn 
00  6 U 4 J = 1 « J c 

8 0 a w R I T F.  ( 8 ♦ 8 0 8 ) I»J»(XlClTtJtJJ)»IJ=l»o) 
bUb  FORMAT  U7t  I 10,F1S.?»5F la.2) 

*s  P I T t.  ( 8 t 8 U 8 ) I t t X 1 T ( J , [ J ) . I J = 1 t O ) 

808  FORMAT  (/'  Total  IN  CONSIST  no. ' t I2»F9.2tSF la. 2//) 

81  C CUNT  J \ 0 a 

v.kiTE  ( 8 ♦ 8 1 ?.  ) 

812  FORMAT  ( 1 HI  t 3SX  t ' I nJuPY  PROFlLt  nY  PASSENGER  C 0 N E I G U r<  A I 1 u N ' / / 1 


00 

6 2 8 

K 

= 1 . a 

IF 

( k . f;  w . 

1 ) 

^1  TE 

( 8 t 8 1 a ) 

IF 

( K • F Li  . 

2) 

R R 1 T E 

l 8 * 8 1 8 J 

IF 

1 K . F 0 . 

5) 

wR  lit 

( 8 t 8 1 8 ] 

IF 

( r 

.El. 

«) 

a1  R I T t 

( 6 t 820  ) 

8 1 a FOkmaT  C////^Ox* ' STANDING  PASSENGFR ' //) 

riltt  FORMA!  (////2Ux.  'FORWARD  FACING  SFATtD  PASSENGER  1 //J 

818  FORMAT  (////20Xf 'SIDEWARD  FACING  SEATED  PASSENGER'//) 

620  FORMAT  ( ////20x  » ' PACKARD  FACING  SEATED  PASSENGER'//) 

8RJTE  (6*822) 

822  FORMAT  ( 

1 3x  t « CONS  I ST  C AR  ' t a 1 X t ' IMJOR Y LtVEL'/6Xt'N0.'*7Xt'N0.'t7X»'uN 

2 N J U R t D ' » 7 x t ‘MINOR ' t 7 x ♦ ' R u D E R A T fc  ' t 7 x t ' S t R 1 0 U S ' ♦ o X » 'CRITICAL' t 5 x t ' 

3 I A L IT] E S ' / ) 

UO  826  1 = 1 *IC 

JC  = NCCI) 

UO  828  J= 1 t J C 

R R 1 T t ( 6 t 8 2 R ) I*J»(xIG(ItJtKtlJ)tTJ=lto) 

- -6  2 4 FORMAT  (17iIlO.F15.2tSr  la. 2) 

626  Continue 


152 


L 

l print  out  - third  phase  (Input  data) 

L 

W w ] T t ( H » c o *'  ) 

9 DU  F On HAT  ( 1 H 1 * 2 9 X * 1 INTERIOR  COUPON ENTS  DATA'/) 

DO  °IN  1 = 1 * I C 
* K 1 T t ( S » 9 0 2 ) ] 

9o?  FORMAT  ( / S X t 'CONSIST  NO.  » .I3//10Xf  ! Jh»  *10X.  »IR»  .15X*  IKC1  » »21X.  IKC^ 
1 • /) 

DU  9o4  jK=l»a 

9 0 9 w R J T fc  ( 6 * 9 U b ) JR*IR(I*JR)*RClfI*JW)*HC2(I*JW) 

9oh  FCjkmaT  (21 1?»?F?2.2) 
a KITE  ( b » 9 u 7 ) 

91)7  F0K--1A7  ( //iX , • PA.-jEL  (JH=1)*  SF.aTdALK  (Jks?),  STAnChIuN  (JP=i)*  FLU 

10  R ( JP  = ii)  • /3*  » • PEPFFC  T AbSOPrE  K (IK=l)t  t L A S I I C RtSPuNSE  (iR=2)»  1 

2 N t L.  A S T I C 9 t S P 0 N S fc  ( 1 * = 3 ) *////) 

W R 7 r F ( h • 9 ( I H ) 

9u*  F 0 K M A 7 ( ////30X  f • c Ap  CRASH  F.  N V I R U N N fc  N T ! / / 1 0 x » ’CONSIST  CAk'tlsX 

l*'vELnr.T7r  Change  ACCELERATION' /12X  t 'no.  inu  . » * l b X . ' ( F 7 /st 

PC ) (FT /SEC  SW)  ' / ) 

DO  910  1=1  »IC 

JC  = N C ( 1 ) 

00  910  J=1*JC 

910  WP17E  (b»912)  I * J» V0( I • J) t A0( I t J) 

91?  FORMAT  (IlUtIlO»F?S.l»Flto.2) 

WRITE  ( b f 9 1 U ) 

914  FORMAT  ( lhl * lbXt 'DISTRIBUTION  FACTORS  bY  CaF  anO  CONSIST'// 

1 3 X . ' COnS 1ST  CAP  STANDING  SE  A T F D FACING  SEATED  FACING1/ 

?b  x ♦ I NO.  ' ♦ 7 X ♦ 'NO.  ' * 1 9 X ♦ 'FORWARD'  * 9X  ♦ 'SIDEwakDI/) 

DO  9 1 b 1=1  t IC 
JC  = MC  ( I ) 

DO  9 1 fc  J=  1 , JC 

91o  aPITE  ( b » 9 1 S ) I f J » (DF ( I f J f K) f K = 1 » 3) 

9 1 F F Ok vat  (I7*llOfFll.2tFlu.2*Fife.?) 

- - - WHITE  (b.920) 

920  FORMAT  (////3*X?  ' SPACING  P A » A * fc T t R S ' / / ) 

DU  92?  1 = 1 * IC 

922  WHITE  (b»92F)  1 ♦ SH ( I ) ♦ SS  ( I ) ♦ b b ( I J ♦ SF  (I) 

929  FORMAT  (bX.  'CONSIST  NO . • ♦ I 3 ♦ 1 0 X ♦ » S H = ' t F b , 2 * SX ♦ ' SS  = ' ♦ F b . 2 ♦ bx » » 5u 
?=  ' ♦F0.2«bX»  »SF  = ' »Fb.?//) 
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<-  o n 


PklUl  U U T 


F ( > U w T h Phi6t 


(OPTIONAL  P K 1 .s  1 OUT  OF  FIXED  DAI  A) 


IF  (IFF  .feU.  U)  UO  Tu  9«e 

W K i T fc  ( p * H *3  0 ) 

PSO  FORMAT  ( 1HJ  »ibx  f 'PARAMETER  VALUES  USED  In  Injury  analy&IS  - FIxE 

1 UAl  A » // ) 

a h ] T fc  ( b » 9 S 2 ) PFintAliiifbnf^HlfPM^tClfCitAntDH^KF^JtL1 1 . A 2 ♦ 


1 w H ♦ P H f 1 ) . P F l *J  } .El 

» C 2 ♦ 

ns. phi?) 

♦ P F l 5 J 

.8  1 

f 

Op  ♦ n 

i.VlFC.bt 

F 0 R * 

AT  ( 

1 1 0*  * 

*PF  ( 1 )r  • .Fa. 2 

♦ 1 0 X 

» ' 1 1 

= ' 

♦ F 5 . 1 v 

1 Ox 

t 

' A3 

= ' » F b . ? ♦ 

21  OX  * 

•00  = ' » F 5 . 2 . 

1 0 X ♦ 

' E>H  1 

r • 

♦ F 5.2/ 

5 H'  x » 

•PF(c)s,»Fa.2 

♦ 1 0 X 

. ' C 1 

r i 

♦ F 5 • 2 * 

1 oX 

f 

' cs 

s ' * F b • 3 * 

a l o x * 

•Ah  - ' t F b « 2 ♦ 

10  x * 

' 8H2 

z • 

* F 5 • ? / 

5 1 0 x . 

'PF(3)= ' .Fu.2 

. 1 o x 

. ' ['1 

r • 

« F 5 , ? * 

1 ox 

f 

' A? 

= ' * F b • 2 » 

r>  1 0 X • 

• h,n  z • » F 5 . 2 . 

1 ox  * 

• PH ( 1 ) 

= • 

* F 5 • 2 / 

7 1 0 X * 

1 P F (9)  = ' t F 4 . 2 

♦ 1 OX 

.'El 

= t 

.Fb.lt 

1 0 x 

t 

' c? 

= » *F5.2* 

81  OX  ♦ 

•I’S  = ' *Fb,2. 

1 ox  * 

' P H ( 2 ) 

= • 

• F b . 2 / 

9 1 0 X * 

'PF(b)='.  Fa. 2 

» 1 0 X 

.'Ft 

— 1 

» F b • 1 ♦ 

1 0 x 

f 

• OP 

= ' * F 5 • 2 * 

n o x . 

'Hi  = ' .Fb,?. 

1 OX  f 

' v I FC 

= ' 

♦Fb.l/box 

t 

' bE 

s ' *Fb.2/) 

1 T fc  ( b » ) 

95P  FUK'IAT  ( ///box  • • SF  VFH]  T Y INDEX  CKI  TERI  A * //5x  f ' THKESnulu  Lfe*vtL**l 
1 ♦ ' hfc AO  • * 21 X ♦ ' F ACE  ' » 2 1 X * ' CHEST  ' / ) 

DU  9bb  I h=  1 t 5 

9bb  *WITE  ((5*958)  18*  (SIC(M.IR)  *«=1  *3) 

958  FORMAT  ( 1 15*3 F 25.0) 

white  ( 6 * 9 o 0 ) SICf«*l) *SIC(A*2)*SIC(«*3) 

9t>0  F 0 h •**  A T ( //20x  * • NECK  INJURY  C R IT F k I A ' . 1 0 X . ' M I nOk  INJURY  THRESHOLD 
1 SIC  = ' .F5.0/50X . ' muDFPATE  INJURY  ThreShuID  * SIC  =»*F5.0/ 
250X*'UPPfcR  L I y I T Cun  I Pul  * SIC  s'*FS.O/) 

998  00  TO  lu 

9 V 9 CONTINUE 
STOP 
End 
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1* 

SUBROUTINE  IMPaCTCJP*SI) 

2* 

COMMON/SUBRT/  I»J*K»VI»AH»wH*IRC2»a)tRCl(2»a),RC2<2 

3* 

IF  (IR(I»K)-2)  1*2*3 

a* 

C 

PERFECT  ABSORBER 

5* 

1 

SI  « ( (RCi  (I*JR)*ah/wh)**i ,5)*Vl/32.2 

6* 

GO  TO  a 

7* 

C 

ELASTIC  PESPONSE 

8* 

2 

SI  = 1 ,a2*( (RCi (If JR)*3»6./WH)**.75)*( (VI/32. 2)**2. 

9* 

GO  TO  a 

10* 

C 

INELASTIC  RESPONSE 

11* 

3 

FN  a RC2(I*JR)/(VI*$QPTCRC1(I*JR)**h*,373)) 

12* 

IF  (FN  ,LT . .1  .OR.  FN  ,GT,  5.)  GO  TO  10 

13* 

F * ( 1 0.** ( ( . 3 1 a * ( ( 1 , + ALOGlO (FN) ) **2,52a}», 15) ) ) 

la* 

GO  TO  20 

15* 

10 

IF  (FN  ,LE.  .1)  F s .71 

lfc* 

IF  (FN  . GE  • 5.)  F s FN**1.5 

17* 

20 

SI  s F*( (PCI (I »JR)*366,/fcrt)**,75)*( (Vl/32,2)**2.5) 

18* 

a 

RETURN 

19* 

END 

1* 

SUBROUTINE  CO*PAR(SIA.m. I J) 

2* 

COMmon/SICRIT/  SIC(a»5) 

3* 

IJ  * 0 

a* 

10 

IJ  ■ IJ+1 

5* 

IF  (IJ  ,EQ.  b)  GO  TO  20 

6* 

IF  ( S I A ,GT.  SIC(M*IJ))  G 

7* 

20 

RETURN 

8* 

ENU 
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interior 

COMPONENTS  DATA 

CONSIST  NO.  1 

JR 

IR 

RC1 

RC2 

1 

3 

20.00 

96,00 

2 

3 

98.00 

307.00 

3 

3 

1 14,00 

218.00 

a 

2 

8000.00 

.00 

CONSIST  NO.  2 

JR 

IR 

RC1 

RC2 

1 

3 

200.00 

300,00 

2 

1 

100,00 

,00 

3 

3 

200.00 

20,00 

a 

2 

2200,00 

.00 

PANEL  ( JR*  1 ) * S&AT5ACK  (JRs2), 

STANCHION  ( J R e 3 ) » 

FLOOR  (JRS4I) 

PERFECT  ABSORBER 

C I R*  1 ) * ELASTIC  RESPONSE  (IR*2) 

. INELASTIC  RESPONSE 

CAR 

CRASH  ENVIRONMENT 

CONSIST 

CAR 

VELOCITY  CHANGE 

ACCELERATION 

NO. 

NO. 

(FT/SEC) 

CFT/SEC  SQ) 

1 

1 

25.2 

5.32 

1 

2 

23. a 

4.37 

1 

3 

22.9 

4.28 

1 

a 

22.8 

4.09 

1 

5 

22,5 

4,08 

1 

6 

22,2 

3.99 

1 

7 

22.1 

3.98 

1 

8 

22.0 

3.97 

2 

1 

41 .3 

7,82 

2 

2 

35.2 

6.41 

2 

3 

34.0 

6.32 

2 

a 

33.7 

6,28 

2 

5 

33.5 

6,24 

2 

6 

32.8 

6.20 
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PARAMETFk  VALUFS  USfcU  IN  INJURY  ANALYSIS  - F I X t D DATA 
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